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Abstract 

The  objective  of  this  study  is  to  investigate  the  effects  of  time-varying,  compressible  flow  on 
rotorcraft  dynamic  stall.  The  problem  was  studied  in  an  unsteady  transonic  wind  tunnel  and 
partitioned  into  three  segments  to  facilitate  understanding:  pitch  oscillations  in  a  steady 
freestream,  investigating  unsteady  effects  (reduced  frequency)  and  compressibility  effects  (Mach 
number);  fixed  pitch  in  a  modulated  freestream,  studying  various  pre-  and  post-stall  angles  with 
fluctuating  Mach  number;  and  oscillations  of  both  angle  of  attack  and  Mach  number,  studying 
the  coupled  effects  and  the  importance  of  phasing  between  the  two.  The  impact  of 
compressibility  on  the  traditional  dynamic  stall  process  was  evaluated,  with  higher  Mach 
numbers  (and  corresponding  increased  compressibility  effects)  changing  the  character  of 
dynamic  stall  on  the  SSC-A09  from  leading-edge  to  a  trailing-edge  phenomenon.  Mach 
oscillations  on  a  fixed-angle  airfoil  were  found  to  have  significant  hysteresis  at  pre-stall  angles, 
and  substantial  unsteadiness  (periodic  pressure  fluctuations)  in  the  post-stall  regime.  Coupled 
pitch  and  Mach  oscillations  showed  that  periods  of  flow  acceleration  have  a  stabilizing  effect  on 
the  dynamic  stall  event. 
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Nomenclature 

A  test  section  area 

A*  throat  area 

c  airfoil  chord 

Cl  lift  coefficient 

CL  lift  curve  slope,  dCJda 

Cm  moment  coefficient 

CMa  moment  curve  slope,  aCMjoa 

Cp  pressure  coefficient 

Cp,cr  critical  pressure  coefficient 
/  frequency,  Hz 

k  reduced  frequency,  coc/lU 

M  Mach  number 

p  pressure 

r  rotor  radial  location 

R  rotor  radius 

Re  Reynolds  number,  Uc/v 

t  time 

T  period 

U  velocity 

x  coordinate  in  the  chordwise  direction 
x'  choke  vane  coordinate 

y  coordinate  perpendicular  to  the  chord 
v'  choke  vane  coordinate 

z  coordinate  in  the  spanwise  direction 
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Greek 


a  angle  of  attack 

ass  static  stall  angle  of  attack 

y  ratio  of  specific  heats 

(p  phase  angle 

Ocv  angular  position  of  the  control  vanes 
p  density 

co  pitching  frequency  (1 /sec) 

C,  distance  along  the  boundary  layer  edge 
Subscripts 

e  edge  of  boundary  layer 

oo  freestream  quantity 
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1.  Introduction  and  Background 

Dynamic  stall  (DS)  is  a  performance-limiting  phenomenon  experienced  by  rotorcraft  in 
forward  flight,  maneuver,  and  at  high  altitude.  The  dynamic  stall  event  typically  occurs  on  the 
retreating  blade  (third  and  fourth  quadrants)  when  the  rotor  blade  is  rapidly  pitched  to  high  angle 
of  attack.  The  lift  stall  and  moment  stall  due  to  the  shed  vortex  can  produce  significant  variations 
in  pitching  moment.  The  transient,  very  high  pitch  link  loads  resulting  from  dynamic  stall  force 
design  choices  that  add  to  the  weight  of  the  vehicle,  and  limit  the  operational  envelope  of  the 
rotorcraft.  Furthermore,  certain  dynamic  stall  phenomena  induce  significant  vibratory  loads  that 
preclude  operation  under  certain  flight  conditions.  These  physical  phenomena  place  severe 
limitations  on  the  maximum  forward  flight  speed  (advance  ratio)  of  rotorcraft,  which  directly 
impacts  Anny  rotorcraft  operations.  In  many  operational  situations,  it  is  very  important  to 
minimize  the  time  required  for  a  helicopter  to  transit  to  theater  operations.  The  maximum 
forward  flight  speed,  which  is  directly  limited  by  compressible  dynamic  stall,  must  be  increased 
in  order  to  improve  deployment  response  times.  Thus,  the  dynamic  stall  phenomenon  directly 
impacts  Anny  aviation  operations. 

A  full,  detailed  understanding  of  the  aerodynamic  phenomena  inherent  to  compressible 
dynamic  stall  can  enable  control  strategies  or  blade  geometry  optimization  that  may  alleviate  the 
very  high  pitch  link  loads  and  expand  the  operating  envelope.  Dynamic  stall  is  a  well-studied 
problem;  however,  most  of  the  historical  investigations  have  not  accounted  for  the  effects  of 
compressibility.  In  particular,  there  is  a  lack  of  knowledge  about  the  effects  of  unsteady 
compressibility  effects  on  dynamic  stall  that  is  typical  of  rotorcraft  operations  in  maneuver  or 
high-speed  forward  flight.  In  order  to  develop  a  full,  accurate  understanding  of  this  dynamic  stall 
phenomenon,  the  problem  must  be  investigated  with  advanced  instrumentation  in  conditions  that 
adequately  model  time-varying  compressibility  effects.  The  purpose  of  this  investigation  is  to 
experimentally  model  the  full-unsteady  compressible  flow  situation  and  perfonn  a  detailed  study 
of  the  complex  compressibility  effects  on  boundary  layer  character  and  DS  vortex  development 
at  matched  reduced  frequency,  Mach  number,  and  Reynolds  number. 

Relevant  studies  of  dynamic  stall  stretch  over  four  decades,  such  that  the  basic  phenomena  of 
dynamic  stall  under  incompressible,  static  freestream  conditions  are  thoroughly  documented 
(McCroskey  et  al.,  1976;  Lorber  and  Carta,  1987;  Lorber,  Carta,  and  Covino,  1992;  Carr,  1988). 
Contemporary  investigations  of  the  intricate  details  of  dynamic  stall  are  also  currently  underway 
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(Geissler  et  al.,  2005;  Sahoo  et  al.,  2008;  Martin  et  al.,  2008;  DiOttavio  et  al.,  2008).  Despite  the 
extensive,  rich  literature  base  in  the  field,  rarely  are  Mach  and  Reynolds  number  varied 
independently  to  study  their  influence  on  DS,  and  compressibility  effects  are  not  well  understood 
(Carr  and  Chandrasekhara,  1996).  Due  to  large  amplitude  pitching  motions,  compressibility 
effects  can  be  important  for  onset  Mach  numbers  as  low  as  0.2.  At  higher  Mach  numbers,  shocks 
can  produce  premature  laminar  separation  during  the  pitch-up  motion  -  thus  interrupting  the 
vortex  shedding  that  is  characteristic  of  deep  dynamic  stall  (Ericsson  and  Reding,  1984).  The 
state  of  the  boundary  layer  and  its  susceptibility  to  separation  are  in  turn  functions  of  Reynolds 
number.  Thus,  the  two  flow  parameters  are  strongly  coupled.  Typical  onset  Mach  and  Reynolds 
number  ranges  are  0.2-0. 5  and  2-6  million,  respectively,  for  retreating  blade  stall.  Mach  numbers 
for  advancing  blade  stall  can  exceed  0.8  (with  correspondingly  higher  Re  values),  as  noted  by 
Bousman  (1998). 

In  order  to  study  the  effects  of  compressibility  on  dynamic  stall,  Carr  and  Chandrasekhara 
(1992)  developed  the  Compressible  Dynamic  Stall  Facility  at  NASA  Ames  Research  Center. 
Their  facility  featured  excellent  optical  access,  freestream  Mach  numbers  up  to  0.5,  and  airfoil 
oscillation  frequencies  up  to  100  Hz,  allowing  reduced  frequencies  up  to  0.15  to  be  tested  in 
compressible  flows.  Their  studies  have  included  determinations  of  the  DS  vortex  convective 
speed  in  compressible  flows  (Chandrasekhara  and  Carr,  1990),  development  of  a  point 
diffraction  interferometry  system  for  quantitative  visualization  of  DS  vortex  behavior  (Carr  et 
al.,  1994),  a  detailed  investigation  of  the  mechanisms  that  drive  compressible  dynamic  stall 
(Chandrasekhara  et  al.,  1998),  multi-element  airfoils  for  DS  alleviation  (Carr  et  al.,  2001),  and 
heat-flux  gauge  investigations  that  captured  the  shock  imprint  near  the  leading  edge 
(Chandrasekhara  and  Wilder,  2003).  One  outcome  of  the  seminal  investigations  into 
compressibility  effects  on  dynamic  stall  is  the  development  of  the  dynamic  drooping  leading 
edge  (Chandrasekhara  et  al.,  2001;  Martin  et  al.,  2003;  Chandrasekhara  et  al.,  2004;  and 
Chandrasekhara,  2007). 

Dynamic  stall  induces  large  amplitude  mechanical  oscillations  of  the  rotor  accompanied  by 
large  pitch-link  loads,  due  to  the  strong  negative  moment  induced  when  the  stall  vortex  travels 
along  the  blade  chord  from  leading  to  trailing  edge.  For  the  retreating  blade  dynamic  stall,  the 
reduced  frequency  of  the  first  harmonic  (1/rev)  blade  oscillation  (k  =  0.5  coc/U)  ranges  from  0.05 
to  0.2.  The  majority  of  helicopter  rotor  dynamic  stall  studies  reported  in  the  literature  have  been 
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conducted  in  a  constant  velocity  wind  tunnel  with  either  a  2D  airfoil  pitching  (or  plunging)  at  the 
appropriate  reduced  frequency.  In  reality,  since  the  rotor  stall  occurs  during  the  retreating  blade 
motion  (at  advance  ratios  that  can  approach  0.3  to  0.4),  the  component  of  the  vehicle  advance 
velocity  “seen”  by  the  retreating  airfoil  also  varies  during  the  dynamic  stall  event.  Thus,  in  a 
rotor-relative  frame,  both  the  approach  flow  angle  (due  to  pitch/plunge)  and  magnitude  are  time- 
varying.  This  time-varying  relative  velocity  has  a  direct  influence  on  the  severity  of  the  local 
pressure  gradient  as  shown  by  Ericsson  (1985).  Ericsson  applied  the  Navier-Stokes  equation  at 
the  boundary  layer  edge  to  obtain  an  expression  for  the  streamwise  pressure  gradient,  dpJdQ. 


J_8p±  =  dU^  +  u  dU^ 

Pe  d£  dt  e  dc 


(1) 


From  this  is  it  clear  that  the  streamwise  pressure  gradient  that  is  imposed  on  the  boundary  layer 
is  a  function  of  both  spatial  and  temporal  gradients  in  velocity.  Namely,  the  adversity  of  the 
pressure  gradient,  dpJdC,,  will  be  decreased  by  an  accelerating  external  velocity,  dUJdt  >  0. 
Since  the  pressure  gradient  is  a  critical  parameter  influencing  boundary  layer  separation,  dUJdt 
will  tend  to  either  promote  or  suppress  separation.  This  led  Ericsson  to  conclude  that,  “The  firee- 
flight/wind  tunnel  equivalence  derived  for  unsteady  inviscid  flow  in  general,  cannot  be  applied  to 
dynamic  stall  analysis  because  of  the  unsteady  viscous  boundary  condition  at  the  airfoil  surface, 
the  so-called  moving  wall  effect.”  (Ericsson,  1985).  Thus,  it  is  evident  that  time-varying 
freestream  velocity  will  have  a  direct  impact  on  the  adversity  of  the  instantaneous  pressure 
gradient,  and  thus,  the  formation  of  the  dynamic  stall  process. 

Several  researchers  have  investigated  the  combination  of  airfoil  pitching  and  relative  velocity 
oscillations  (Pierce  et  al.,  1978;  Favier  et  al.,  1988)  and  concluded  that  a  strong  coupling  exists 
between  the  two  modes  of  oscillation.  Of  note,  all  research  of  this  type  has  been  conducted  at 
incompressible  flow  speeds  below  Mach  0.1.  In  contrast,  the  effort  presented  here  studies  the 
effects  of  unsteady  Mach  oscillations  on  the  dynamic  stall  process,  where  the  unsteady 
compressibility  effects  are  largely  unknown. 

Fernie  and  Babinsky  (2003,  2004)  conducted  a  very  interesting  investigation  into  the  effects 
of  unsteady  freestream  oscillations  on  the  dynamic  shock  structure  on  a  NACA  0012  airfoil. 
When  operated  at  transonic  speeds  and  moderately  high  frequencies  (of  the  order  of  50  Hz),  they 
found  that  there  were  differences  in  the  shock  strength  and  speed,  depending  on  whether  the 
freestream  Mach  number  was  increasing  or  decreasing.  They  found  that  a  decelerating 
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freestream  caused  a  stronger  shock  to  move  forward  on  the  airfoil  surface  more  rapidly  than  an 
accelerating  freestream  over  the  same  change  in  Mach  number.  Fernie  and  Babinsky’s  work 
shows  that  compressibility  effects  should  be  significant,  yet  they  stopped  short  of  perfonning  any 
dynamic  stall  investigations  (their  work  was  limited  to  a  fixed  airfoil). 

The  work  presented  here  is  unique  in  that  it  properly  models  the  effects  of  time-varying 
Mach  number  on  dynamic  stall  phenomena  at  suitable  reduced  frequencies  and  advance  ratios 
typical  of  high-speed  forward  flight  of  rotorcraft.  The  objectives  of  this  experimental  research 
include  establishing  the  significance  of  and  understanding  the  impact  of  time-varying 
compressibility  on  compressible  dynamic  stall.  There  are  three  unsteady  scenarios  which  are 
simulated  experimentally.  The  first  scenario  is  that  of  the  oscillating  airfoil  in  a  steady  airstream. 
The  dynamic  stall  encountered  in  this  scenario  constitutes  the  conventional  simulation  of  the 
dynamic  stall  process  that  has  contributed  to  our  understanding  of  the  key  physical  phenomena. 
The  second  scenario  involves  evaluation  of  how  hysteresis  and  stall  are  impacted  when  a  fixed- 
pitch  airfoil  is  subjected  to  a  time-varying  freestream  Mach  number  condition.  Finally,  the  airfoil 
model  oscillations  are  synchronized  with  the  freestream  velocity  oscillations  for  coupled  pitch 
and  Mach  oscillations.  Moreover,  measurements  are  made  while  varying  the  phasing  relationship 
between  these  synchronized  oscillations.  Through  this  progression,  the  impact  of  time -varying 
compressibility  on  the  compressible  dynamic  stall  process  is  elucidated. 

2.  Experimental  Configuration 

2.1  Transonic  Wind  Tunnel 

The  dynamic  stall  problem  was  investigated  in  the  6”  x  22”  transonic  wind  tunnel  located  in 
the  Aeronautical  and  Astronautical  Research  Laboratories  (AARL)  at  The  Ohio  State  University. 
The  facility  has  been  used  over  the  last  three  decades  by  every  major  aircraft  manufacturer  for 
airfoil  studies  (Lee  et  ah,  1978).  The  tunnel  has  also  been  used  for  dynamic  stall  investigations 
on  airfoil  sections  such  as  the  NACA  0015  (Gregorek  et  ah,  1989).  A  layout  of  the  tunnel  can  be 
seen  in  Figure  1.  Inlet  flow  to  this  open-circuit  wind  tunnel  is  delivered  by  20cm  (8in.)  supply 
line  from  two  42. 5nr  (1500  ft  )  air  storage  tanks  pressurized  up  to  17MPa  (2500  psi)  with  in-line 
air  dryers  to  maintain  gas  purity.  The  high  pressure  air  flow  is  controlled  by  two  valves.  The  first 
is  a  control  valve  which  sets  the  total  pressure  and  Reynolds  number.  The  second  valve  is  a  fast 
acting  valve  used  to  start  and  stop  the  flow.  Valving  and  various  multifunction  data  acquisition 


J.  W.  Gregory,  W91 1NF1 1 10503  Final  Report,  p.  13 


tasks  are  orchestrated  using  the  6”  x  22”  control  computer.  The  maximum  operating  pressure  of 
the  wind  tunnel  is  350kPa  (50psia).  The  settling  chamber  is  equipped  with  a  perforated  plate, 
honeycomb  section  and  then  eight  screens  (60-mesh)  to  condition  the  flow  and  lower  the  test- 
section  turbulence  intensity  to  <  0.5%  of  the  freestream  velocity  under  steady  flow  conditions.  A 
subsonic  nozzle  with  a  contraction  ratio  of  15:1  provides  excellent  flow  uniformity  in  the  6”  x 
22”  test  section  which  is  1.1m  long.  The  56cm  (22in.)  high,  solid  sidewalls  hold  the  airfoil.  The 
15cm  (6in.)  spanwise  floor  and  ceiling  walls  are  perforated  with  3.2mm  straight  holes  yielding 
an  effective  porosity  of  6  percent  (Petrie  and  Davis,  1978).  These  isolation  cavities  or  plenums 
are  open  to  the  flow  only  downstream  of  the  test  section  and  aid  in  producing  a  high  quality  flow 
in  the  test  section  by  reducing  Mach  wave  reflections  in  transonic  flow. 

High  Pressure  Air 
(42.5  mJ  at  15.5MPa) 


Figure  1:  Schematic  of  the  Aeronautical  &  Astronautical  Research  Lab  6”  x  22”  transonic 

wind  tunnel. 

Mach  number  is  controlled  independently  of  Reynolds  number  by  adjusting  the  throat  area 
downstream  of  the  test  section.  The  Mach  number  is  uniquely  set  by  the  throat  area,  independent 
of  stagnation  pressure,  as  long  as  choked  flow  is  maintained.  The  two  parameters  (Re  and  M)  can 
be  independently  varied  over  a  considerable  range,  as  shown  in  Figure  2. 
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Figure  2:  Range  of  operating  conditions  for  the  6”  x  22”  transonic  wind  tunnel  at  287K 
(56°F)  and  a  max  operating  pressure  of  345kPa  (50psia). 

The  6”  x  22”  is  equipped  to  operate  in  several  dynamic  modes  in  order  to  simulate 
compressible  dynamic  stall.  Its  accessories  include  an  airfoil  oscillation  mechanism  which  is 
based  on  a  5-hp  motor  that  drives  a  face  cam  and  mechanical  link  assembly,  as  shown  in  Figure 
3  through  Figure  5.  In  its  current  configuration,  the  oscillation  mechanism  has  been  reliably 
operated  at  frequencies  up  to  17  Hz  over  a  wide  range  of  oscillation  amplitudes.  Also,  the  6”  x 
22”  transonic  wind  tunnel  facility  was  recently  modified  by  the  installation  of  rotating  vanes  at 
the  throat,  to  allow  for  variable  choke  area  and  Mach  oscillations  at  rates  up  to  17  Hz  in  the  test 
section  (Gompertz  et  ai,  2011).  An  overview  of  the  airfoil  pitch  oscillation  and  Mach  oscillation 
devices  is  presented  in  Figure  3  and  they  are  discussed  in  more  detail  in  subsequent  sections. 
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Figure  3:  Overview  of  Mach  oscillation  and  pitch  oscillation  drive  mechanisms. 

For  an  airfoil  with  a  152.4mm  (6”)  chord  length,  the  Mach  number  range  (0.2- 1.0),  Reynolds 
number  range  (1-8  million),  and  unsteady  reduced  frequency  capabilities  of  this  facility  are 
ideally  suited  for  studying  the  rotor  dynamic  stall  phenomenon  for  both  the  retreating  and 
advancing  blade. 

2.2  Airfoil  Pitch  Oscillation  Assembly 

An  airfoil  oscillation  mechanism  is  installed  in  the  6”  x  22”  transonic  wind  tunnel  as  seen  in 
Figure  3  through  Figure  5.  The  airfoil  oscillates  about  the  quarter  chord  position  by  a  5-hp 
variable-frequency  A/C  motor.  The  drive  chain  is  configured  such  that  a  belt  rotates  an  eccentric 
disk  and  a  momentum  fly  wheel  to  minimize  the  unsteady  aerodynamic  loads  on  the 
transmission.  A  lever  arm  is  connected  to  the  eccentric  disk  by  a  pin,  a  connection  bar  then 
oscillates  the  airfoil  which  has  bearings  isolating  it  from  the  tunnel  walls.  For  this  experimental 
study,  the  pitch  oscillation  amplitude  (fixed  by  the  groove  in  the  eccentric  disk)  can  be  set  to 
either  ±5°  or  ±10°.  The  mean  angle  of  attack  is  set  by  the  connection  bar  length.  The  angle  of 
attack  is  quickly  adjusted  from  0°  to  15°  in  increments  of  5°  using  interchangeable  links 
fastening  the  top  of  the  connection  bar  to  the  airfoil  axis  of  rotation  (Figure  5(b)).  The  oscillating 
test  section  has  a  frequency  range  of  1-17  Hz.  With  this  range,  it  is  possible  to  achieve  an 
unsteady  reduced  frequency  of  approximately  0.05  at  Mach  numbers  up  to  0.45.  A  shaft- 
mounted  optical  encoder  with  quadrature  outputs  at  2500  counts  per  revolution  is  coupled  to  the 
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shaft  of  the  eccentric  disk.  This  signal  is  sampled  at  50  kHz  and  used  to  detennine  the  angular 
position  of  the  airfoil  to  within  ±0.05°. 


Pulley  diameter  ratio  sets  frequency  range 
Figure  4:  Schematic  of  the  airfoil  oscillation  device. 


(»)  (b) 

Figure  5:  Detail  of  (a)  eccentric  disk  and  drive  mechanism,  and  (b)  thrust  bearing  for  pitch 

oscillation. 
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2.3  Setting  Test-Section  Mach  Number 

In  the  conventional,  time-invariant  configuration,  test-section  Mach  number  is  uniquely 
controlled  independently  of  Reynolds  number  by  adjusting  the  throat  area  (A  )  downstream  of 
the  test  section.  Since  the  blowdown  tunnel  is  choked  at  the  throat,  the  test  section  Mach  number 
is  uniquely  set  by  the  area  ratio  independent  of  the  stagnation  pressure  according  to  the  area- 
Mach  number  relation  expressed  in  Eq.  (2).  The  throat  area  is  set  by  installing  choke  bars  of 
varied  diameter  as  indicated  in  Figure  1 . 
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2.4  Freestream  Mach  Oscillation  Assembly 

The  transonic  tunnel  has  been  modified  to  allow  freestream  velocity  modulation  up  to  17Hz 
using  a  variable-area  throat  downstream  of  the  test  section  (refer  to  Figure  3).  The  modification 
uses  rotating  elongated  vanes  to  control  the  throat  cross-sectional  area  as  a  function  of  time.  A 
drive  mechanism  rotates  four  vanes  of  the  profile  displayed  in  Figure  6.  As  the  vanes  rotate,  the 

open  cross-sectional  area  is  varied  between  maximum  and  minimum  values,  as  illustrated  by  the 

•  •  •  • 
time-history  of  A  also  included  in  Figure  6.  The  cross-sectional  area  of  the  vanes  was  specified 

such  that  the  vane  rotation  produces  a  sinusoidal  Mach  number  oscillation  (0.4  ±  0.08)  in  the  test 

section.  The  amplitude  of  Mach  number  oscillation  is  set  by  the  dimensions  of  the  major  and 

minor  axes  of  the  rotating  vanes,  and  the  specific  Mach  waveform  dictated  by  the  geometry  of 

the  vanes. 

Design  of  the  choke  vanes  for  a  specific  Mach  waveform  (say,  M  =  0.4  +  0.08  sin  9  in  this 
case)  is  not  straightforward  because  there  is  no  detenninistic  mapping  between  a  required  A* 
time  history  and  a  parameterized  vane  geometry.  Thus,  the  necessary  parameterized  vane  profile 
for  a  given  A  history  must  be  detennined  iteratively.  In  order  to  make  this  iterative  process 
tractable,  a  superellipse  profile  is  used  to  define  the  vane  geometry.  The  generic  definition  of  a 
superellipse  is 
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where  m  and  n  are  both  greater  than  0.  Parametric  equations  for  the  vane  geometry  (x'(0),  y\B )) 
in  Cartesian  coordinates  are  given  by 
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(4) 


,  ,2/ 

x(6,j  =  |cos6,|  m  -flsgnlcos^j 
2/ 

sin#  n  -bsgn^sin#) 

the  major  (a)  and  minor  ( b )  axes  of  the  profile  determine  the  maximum  and  minimum  values  of 
the  Mach  waveform.  Each  of  the  four  vanes  used  in  this  work  have  a  =  2.07  cm  (0.82  in.)  and  b 
=  3.25  cm  (1.28  in.).  A  nonlinear  regression  was  performed  while  allowing  the  regression 
variables  m  and  n  to  vary  to  detennine  the  geometry  which  produces  a  prescribed  sinusoidal  test- 
section  Mach  number  when  rotated  at  constant  angular  velocity.  The  parameters  which  define 
this  choke  vane  geometry  are  m  =  1.9833  and  n  =  1.5402. 

The  geometrical  configuration  of  the  oscillating  vanes  may  be  adjusted  between  runs  in  a 
relatively  short  period  of  time  by  interchanging  between  several  sets  of  pre-fabricated  vane 
assemblies  to  produce  varied  Mach  number  waveforms.  The  vane  assembly  is  driven  by  a  4  kW 
A/C  servo  motor,  which  is  slaved  to  the  pitch  oscillation  motor.  The  drive  chain  is  configured 
such  that  adjacent  vanes  are  counter-rotating,  to  minimize  the  unsteady  aerodynamic  loads  on  the 
transmission.  A  shaft-mount  optical  encoder  is  coupled  to  the  shaft  of  the  upper-most  choke  vane 
to  index  each  full  revolution  (1/rev)  and  the  angular  orientation  of  the  choke  vanes  (500/rev). 
The  signal  is  used  to  determine  the  oscillation  frequency  and  the  instantaneous  angular  position 
(0cv). 


Throat  Area  \s.  Choke  Vane  Angular  Position 
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Figure  6:  Choke  vane  geometry  and  resulting  physical  throat  area  variation. 
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2.5  Freestream  Mach  Number  Measurement 

It  is  well-known  that  modeling  transonic  airfoil  stall  in  a  wind  tunnel  can  be  challenging  due 
to  confinement  of  streamlines  that  would  deform  differently  in  a  real  freestream  condition 
(Goethert,  2007).  Goethert  suggests  suitable  freestream  velocity  corrections  that  would  be 
appropriate  for  this  work.  However,  it  is  important  to  select  a  suitable  measurement  station  for 
the  freestream  Mach  number.  The  location  should  ensure  that  the  airfoil  flow  is  not  perturbed  by 
the  presence  of  the  probe,  and  the  potential  field  distortion  that  results  from  the  presence  of  the 
airfoil  should  be  negligible.  The  potential  field  distortion  is  highest  when  the  airfoil  is  at  high 
angle  of  attack,  which  is  the  condition  of  interest  for  this  type  of  study.  Furthermore,  the  time- 
varying  freestream  velocity  introduces  another  very  important  factor  in  determining  the 
measurement  station.  In  the  time-invariant  case,  the  freestream  static  pressure  is  measured  within 
the  perforated  plenums  (Figure  1)  located  above  and  below  the  test  section.  Ordinarily,  these 
plenums  serve  to  better  simulate  a  freestream  condition,  and  dampen  shock  waves  that  would 
otherwise  reflect  and  possibly  contaminate  the  test-section.  However,  in  the  case  of  unsteady 
flow,  the  dynamics  in  these  plenums  distort  the  measured  freestream  velocity  wavefonn. 

Two-dimensional  computational  simulations  were  performed  with  Fluent  12.0  to  determine 
the  extent  of  the  upstream  influence  of  an  airfoil  installed  in  the  6”  x  22”,  and  to  aid  in  selecting 
a  suitable  location  to  place  a  probe  for  freestream  Mach  number  measurement.  These  simulations 
covered  a  range  of  freestream  Mach  numbers  from  0.3  to  0.65.  The  tunnel  was  modeled  as  a  long 
rectangular  duct  with  solid  wall  boundaries.  The  simulations  included  viscous  effects  to  capture 
first  order  shock-boundary  layer  effects,  and  the  Spalart-Allmaras  turbulence  model  was  used. 
The  simulations  were  validated  by  comparing  predicted  Cp  distributions  with  experimental  data 
for  a  couple  of  airfoils  (but  not  the  SSC-A09  used  in  this  study).  Also  sidewall  static  pressure 
distributions  were  compared  to  ensure  validity  of  the  prediction  tool  for  this  purpose.  Figure  7 
shows  the  domain  and  the  location  of  data  traces  that  are  presented  in  Figure  8.  Observing  Figure 
8  for  a  sample  case  where  the  inlet  Mach  number  is  approximately  0.3,  it  is  evident  that  placing 
the  probe  0.3  to  0.5m  upstream  of  the  airfoil,  and  0.2m  off  of  the  mid-height  will  result  in  a 
measured  Mach  number  with  negligible  distortion  due  to  potential  effects.  Similar  conclusions 
are  drawn  for  other  values  of  freestream  Mach  number  and  airfoil  angle  of  attack. 
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Transverse  Traces 


A  Pitot-static  probe  is  mounted  approximately  0.5m  (3.3  chord  lengths)  upstream  of  the 
airfoil  leading  edge  and  0.2m  (1.3  chord  lengths)  above  the  airfoil.  Each  sensor  orifice  is 
connected  to  a  miniature  pressure  transducer  via  a  short  length  of  plastic  tubing.  Pneumatic 
corrections  are  applied  in  the  same  manner  outlined  in  a  later  section  discussing  dynamic 
response  considerations  for  unsteady  pressure  measurements.  The  Pitot-static  probe  accurately 
resolves  pressure  fluctuations  at  frequencies  up  to  500Hz  which  is  sufficient  to  resolve  the 
expected  the  frequency  content.  The  ratio  of  static  pressure  to  the  stagnation  pressure  is  used  to 
calculate  the  test  section  Mach  number  via  the  isentropic  pressure  relationship.  The  calibration 
uncertainty  of  measured  Mach  number  is  ±0.005. 
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Transverse  Mach  Number  Traces 


Figure  8:  Transverse  traces  of  predicted  Mach  number  at  several  locations  upstream  of  test 

article. 

2.6  Stagnation  Temperature 

The  flow  temperature  decreases  at  a  rate  of  approximately  0.5K/s  (l°F/s)  as  air  is  discharged 
from  the  constant-volume  high  pressure  tanks.  Consequently,  the  speed  of  sound  can  decrease  by 
as  much  as  9  m/s  during  a  30-second  run.  A  Type  J  thermocouple  is  used  to  measure  the 
stagnation  temperature  which  is  recorded  at  a  rate  of  10  Hz  to  the  6”  x  22”  control  computer. 
Static  temperature  in  the  test  section  is  detennined  via  the  isentropic  temperature  relation  using 
the  freestream  Mach  number  (measured  with  the  Pitot-static  probe). 

2. 7  Airfoil  Surface  Pressure  Data 

Airfoil  surface  pressure  taps  are  connected  to  an  ESP  Pressure  Scanner  System  to  obtain  the 
static  pressure  distribution  along  the  airfoil  surface.  The  pressure-coefficient  (Cp)  is  calculated  at 
each  tap  location  as  defined  in  Eqn.  (5), 
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The  trapezoidal  rule  is  used  to  integrate  the  Cp  curve,  and  calculate  the  sectional  lift  (C£)  and 
moment  (Cm)  coefficients.  It  is  noted  that  C/  and  Cm  measurements  inherently  omit  contributions 
of  skin  friction  forces. 

Two  ESP  32HD  Pressure  scanners  are  mounted  just  outside  the  wind  tunnel  and  the  pressure 
ports  are  connected  to  airfoil  pressures  via  durable  plastic  tubing  of  lmm-diameter  (0.040in.). 
The  ESP  scanners  are  miniature  electronic  differential  pressure  measurement  units  consisting  of 
an  array  of  32  silicon  piezoresistive  pressure  sensors  (±210kPa  differential),  one  for  each  port. 
The  electrical  outputs  of  the  sensors  are  electronically  multiplexed  through  a  single  onboard 
instrumentation  amplifier  at  switching  rates  up  to  70,000Hz  using  binary  addressing.  The 
multiplexed  amplified  analog  output  is  connected  to  a  DTC  Initium  scanner  interface  and 
streamed  via  Ethernet  to  a  dedicated  hard  disk.  In  order  to  ensure  data  streaming  at  a  stable 
sample  interval,  the  Initium  is  triggered  with  an  external  TTL  pulse  train  which  is  generated  by 
the  6”  x  22”  control  computer.  In  this  hardware  trigger  mode,  the  maximum  scanning  rate  is 
1,200  samples/sec  per  sensor.  Also,  triggering  assures  accurate  temporal  correlation  with  the 
various  other  analog  voltage  signals  being  recorded  by  the  control  computer.  Aside  from  the  high 
data  rates,  the  ESP  Scanner  System  features  individual  sensor  thennal  compensation  to  minimize 
zero  and  span  shifts  with  temperature  which  are  inherent  in  silicon  pressure  sensors. 

The  ESP  scanners  also  incorporate  a  two-position  calibration  manifold  actuated  by 
momentary  pulses  of  control  pressures.  When  placed  in  the  calibrate  position,  all  sensors  are 
connected  to  a  common  calibration  pressure  port.  To  ensure  the  highest  accuracy,  an  automated 
five-point  calibration  is  performed  periodically  by  applying  a  sequence  of  accurately  measured 
pressures  through  this  port.  The  transducers  are  accurate  to  within  ±0.03%  or  better  of  the  full 
scale  (FS)  pressure  range. 

2.8  Dynamic  Response  Considerations 

Combined  with  the  effect  of  tubing  diameter  and  length,  the  small  ESP  sensor  volume  (with 
high  natural  frequencies)  is  key  to  minimizing  the  pneumatic  distortion  associated  with  remote 
pressure  sensing  schemes.  A  minimal  length  of  tubing  (approximately  20cm  long)  was  used  to 
uniformly  transmit  all  airfoil  surface  pressures  to  the  ESP  Pressure  Scanners. 
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Attenuation  of  measured  pressure  signals  due  to  viscous  effects  in  the  tubing-sensor  system 
is  considered  before  taking  dynamic  pressure  measurements.  Bergh  and  Tij deman  developed  an 
analytic  model  which  corrects  for  attenuation  and  phase  lag  associated  with  pneumatic  tubing 
(Bergh  and  Tijdeman,  1965).  The  model  characterizes  the  dynamic  response  of  the  tubing  as  a 
function  of  the  geometry  of  the  tubing  (inner  diameter  and  length)  and  transducer  (sensor 
volume),  along  with  the  ambient  conditions.  The  Bergh  and  Tijdeman  model  (referred  to  as  B-T 
in  Figure  9)  is  a  transfer  function  that  ratios  the  pressure  measured  by  the  transducer  with  the 
actual  unsteady  pressure  at  the  surface  pressure  tap  in  the  frequency  domain.  An  acoustic 
excitation  test  bed  was  used  to  generate  sinusoidal  pressure  oscillations  over  a  relevant  frequency 
range.  The  oscillating  pressure  was  measured  with  representative  tubing  connecting  pressure 
ports  in  the  test  bed  to  the  ESP  scanners.  Also,  the  oscillating  pressure  was  measured  with  a 
flush-mounted  piezoelectric  transducer.  The  amplitude  ratio  (AR  =  amplitude  measured  with 
ESP  :  amplitude  measured  with  Kulite)  is  plotted  vs.  excitation  frequency  in  Figure  9  along  with 
the  B-T  model  predictions.  The  empirical  data  and  model  predictions  are  in  excellent  agreement, 
providing  a  suitable  scheme  for  attenuation  and  phase  compensation  of  unsteady  pressure  data. 
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Tubing  Diameter,  d  =  1  mm  (0.040  in.) 


Curves  Parametrized  by  Tube  Length 
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Figure  9:  Comparison  of  measured  ESP  pressure  sensor  dynamic  response  and  Bergh- 
Tijdeman  model  (ref.  Bergh  and  Tijdeman,  1965)  for  representative  lengths  of  pneumatic 

tubing  (diameter  =  1  mm  =  0.040in.). 


2.9  Test  Article 

The  model  used  in  this  experimental  campaign  is  an  airfoil  with  the  Sikorsky  SSC-A09 
profile  which  is  shown  in  Figure  10.  The  surface  coordinates  and  derivatives  are  identical  to 
those  published  by  (Lorber  and  Carta,  1987)  and  they  are  tabulated  in  Appendix  1.  The  leading 
edge  radius  is  0.7%  of  chord.  The  maximum  thickness  of  the  airfoil  is  9%  located  near  38% 
chord.  The  thickness  distribution  is  relatively  flat  between  20%  and  60%  chord  to  promote 
healthy  performance  under  supercritical  flow  conditions. 

The  model  is  fabricated  from  solid  aluminum  using  wire  Electrical  Discharge  Machining.  An 
isometric  view  of  the  airfoil  drawing  is  displayed  in  Figure  10(a).  Both  the  chord  and  the  span 
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are  152.4mm  (6  in.),  yielding  an  aspect  ratio  equal  to  1.  One  end  of  the  model  mates  with  a 
“rectangular  side  block”  which  is  used  to  rigidly  mount  the  airfoil  in  the  wind  tunnel  test  section 
while  the  other  end  terminates  with  a  “tap  output  boss”.  The  model  has  a  25.4mm-deep  pocket 
(with  a  minimum  wall  thickness  of  1.27mm)  carved  into  the  boss-side  wherein  all  the  plastic 
lines  transmitting  model  surface  pressures  are  bundled  and  snaked  through  the  tap-output  boss. 

The  model  is  outfitted  with  53  surface  pressure  taps  to  measure  aerodynamic  forces  and 
moments.  Tap  locations  are  shown  in  Figure  10(b)  and  tabulated  in  Appendix  2.  There  are  30 
taps  on  the  upper  surface,  and  23  on  the  lower  surface  (including  the  leading  edge  tap).  The  tap 
distribution  covers  the  leading  82%  of  chord.  Coverage  on  the  trailing  20%  is  challenging 
because  of  the  limited  workspace  within  the  pocket  design  and  the  airfoil  thickness  in  general.  In 
the  absence  of  a  trailing  edge  tap,  a  ghost  tap  is  created  in  post-processing  whereby  the  trailing 
edge  pressure  is  estimated  to  be  the  average  of  the  downstream-most  taps  on  the  upper  and  lower 
surfaces.  Tap  spacing  is  denser  in  the  leading  30%  of  chord  in  order  to  better  resolve  the  large 
pressure  gradients.  However,  the  spatial  resolution  of  surface  pressure  within  the  leading  2.5% 
chord  is  limited  by  manufacturing  and  design  constraints.  The  internal  spanwise  cavities  which 
transmit  surface  pressures  to  the  edge  of  the  airfoil  are  separated  by  one  hole  diameter  of 
material  to  avoid  breaching  between  adjacent  cavities  during  manufacturing.  The  hole  diameter 
is  limited  by  the  distance  the  drill  must  plunge  to  reach  the  midspan.  An  assessment  of  the 
impact  of  tap  resolution  on  the  integrated  quantities  of  lift  and  moment  was  given  by  Hird  et  al. 
(2014). 

The  taps  are  staggered  across  the  middle  third  of  the  model  span  to  prevent  adjacent  taps 
from  interfering  with  one  another.  All  taps  are  located  on  a  plane  which  makes  a  19°  angle  with 
the  airfoil  midspan  plane.  The  stagger  angle  is  such  that  the  turbulent  wedges  which  might 
originate  from  a  pressure  tap  do  not  cover  other  taps.  Taps  are  plunged  normal  to  the  surface 
with  a  0.3mm-diameter  (0.013in.)  to  a  depth  of  1.27mm  where  they  intercept  the  1.27mm- 
diameter  (0.050in.)  spanwise  pressure  lines.  All  internal  spanwise  pressure  lines  in  this  model 
tenninate  in  the  airfoil  pocket  with  a  stainless  steel  tubulation  with  1 ,27mm  inner  diameter  that 
connect  via  durable  plastic  tubing  to  the  pressure  scanners.  The  exit  orifices  can  be  seen  on  one 
side  in  Figure  10.  Uniform  tap  geometry  was  used  in  an  effort  to  simplify  the  dynamic  response 
characterization  and  compensation  which  is  discussed  later.  Thus,  the  total  length  of  the  tubing 
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measured  from  each  tap  to  the  pressure  transducer  is  20cm  (8in.),  and  there  are  no  discontinuities 
of  tubing  inner  diameter  (1.27mm)  along  that  length. 


Staggered  Pressure  Taps 

Angle  w.r.t.  streamwise  dir.  =  24deg 


(a) 


(b) 

Figure  10:  SSC-A09  Airfoil  used  in  this  work:  (a)  isometric  view  and  (b)  side  view  showing 

tap  distribution. 


2.10  Uncertainty  Estimates 

An  analysis  was  conducted  in  the  manner  outlined  by  Coleman  and  Steele  (1999)  to  estimate 
the  relevant  calibration  uncertainties  at  M=  0.3  with  the  wind  tunnel  in  steady-flow  mode:  Mach 
number,  +/-0.005;  Reynolds  number,  ±5,000,  angle  of  attack,  ±0.05°,  Cp,  ±0.05,  Cl,  ±0.05,  Cm, 
±0.02,  CD,  ±0.001,  and  TM,  ±0.05K. 


3.  Results 

The  objectives  of  this  experimental  research  include  establishing  the  significance  of  and 
understanding  the  impact  of  time-varying  compressibility  on  compressible  dynamic  stall.  There 
are  three  unsteady  scenarios  which  are  simulated  experimentally.  The  first  scenario  is  that  of  the 
oscillating  airfoil  in  a  steady  airstream.  The  dynamic  stall  encountered  in  this  scenario  constitutes 
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the  conventional  simulation  of  the  dynamic  stall  process  that  has  contributed  to  our 
understanding  of  the  key  physical  phenomena.  The  second  scenario  involves  evaluation  of  how 
hysteresis  and  stall  are  impacted  when  a  fixed-pitch  airfoil  is  subjected  to  a  time-varying 
freestream  Mach  number  condition.  Finally,  the  airfoil  model  oscillations  are  synchronized  with 
the  freestream  velocity  oscillations  for  coupled  pitch  and  Mach  oscillations.  Moreover, 
measurements  are  made  while  varying  the  phasing  relationship  between  these  synchronized 
oscillations.  Through  this  progression,  the  impact  of  time-varying  compressibility  on  the 
compressible  dynamic  stall  process  is  elucidated. 

3. 1  Pitch  Oscillations  with  Fixed  Mach  number 

The  first  phase  of  this  work  is  focused  on  the  conventional  understanding  of  dynamic  stall  for 
the  case  of  an  airfoil  model  executing  sinusoidal  oscillations  of  pitch  in  a  steady  airstream. 
Comparisons  are  made  with  historical  data  (Lorber  and  Carta,  1987)  for  M  =  0.2  and  a  ~  10°  - 
10°  cos  cp  to  validate  the  experimental  facility  and  methodology.  Then  the  impact  of 
compressibility  on  dynamic  stall  is  systematically  inspected  forM=  0.2,  0.4,  and  0.6  and  a  ~  10° 
-  10°  cos  <p. 

3.1.1  Introductory  Comments 

Figure  1 1  is  a  composite  figure  which  shows  phase-averaged  measurements  for  a  typical 
dynamic  stall  test  in  the  OSU  6”  x  22”  transonic  wind  tunnel  for  the  case  of  an  oscillating  airfoil 
in  a  steady  airstream.  The  particular  test  conditions  are  M  =  0.2,  a  ~  9°  -  11°  cos  cp,  and  k  = 
0.050.  As  with  all  pitch  oscillation  data  reported  in  the  current  study,  maximum  a  is  set  to  20° 
with  the  pitch  oscillation  assembly  unpowered.  The  figure  includes  curves  of  lift  and  moment 
coefficient  on  the  left-hand-side.  Contours  of  upper  surface  pressure  coefficient  are  plotted  in  the 
middle  with  chordwise  position  along  the  abscissa  and  phase  angle  along  the  ordinate.  The  angle 
of  attack  is  plotted  on  the  right  with  angle  along  the  abscissa  and  phase  angle  along  the  ordinate. 
The  Cp  contours  show  the  development  of  suction  during  the  upstroke  tenninating  with  the 
formation  and  convection  of  the  leading-edge  vortex  (LEV).  The  tongue  of  depressed  Cp 
contours  which  extends  over  the  leading  50%  of  chord  from  140°  <  cp  <  160°  is  the  LEV 
signature.  As  the  LEV  convects  past  the  trailing  edge  of  the  airfoil,  the  Cp  plateaus  for  160°  <  9 
<  225°,  indicative  of  an  open  separation  starting  from  the  leading  edge.  This  persists  while  the 
airfoil  begins  the  downstroke  until  cp  =  225°  when  reattachment  begins  from  the  leading  edge  and 
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progresses  downstream  toward  the  trailing  edge  until  the  flow  is  fully  attached  at  (p  ~  270°. 
Integrating  the  phase-averaged  Cp  distribution  in  this  sequence  produces  the  hysteresis  CL  and 
Cm  loops  presented  on  the  left.  While  suction  development  is  stable  during  the  upstroke  the  lift- 
curve  slope  is  approximately  constant  until  a  ~  18°.  The  elevated  suction  coinciding  with  the 
LEV  formation  can  produce  a  brief  dynamic  lift  enhancement  before  the  vortex  convects  past  the 
airfoil.  Evidence  of  this  lift  enhancement  is  sometimes  more  pronounced  for  different  cases, 
though  not  to  the  expected  level  which  will  be  discussed  later.  The  lift  coefficient  abruptly 
decreases  from  1.5  to  0.8  for  19°  <  a  <  20°  while  the  vortex  passes  over  the  airfoil  surface.  As 
the  downstroke  begins  and  the  airfoil  is  completely  stalled,  Cl  is  at  a  value  of  approximately  0.6 
until  a  =  8.5°  {(p  =  270°)  whereafter  the  lift-curve  slope  matches  the  value  measured  during  the 
upstroke.  The  moment  coefficient  is  of  special  interest  when  analyzing  the  dynamic  stall  process 
because  of  the  aeroelastic  implications.  The  impulsive  aerodynamic  moment  experienced  by  a 
rotor  is  transmitted  through  the  structure  to  the  pitch  linkage  at  the  rotor  hub.  By  virtue  of  the 
coupling  between  the  aerodynamic  and  structural  dynamics,  this  cyclic  forcing  can  cause  flutter 
and/or  high-cycle  fatigue  of  the  pitch  linkage.  Thus,  it  is  important  that  the  measurement 
adequately  resolve  the  pitching  moment  behavior.  During  the  upstroke,  for  -2°  <  a  <  18°,  CM 
does  not  vary  suddenly  although  it  has  as  slightly  negative  slope  which  will  be  discussed  in  more 
detail.  When  the  LEV  is  initialized,  and  before  it  convects  past  the  quarter  chord,  the  pressure 
near  the  LE  is  lower,  and  the  moment  about  the  quarter  chord  is  slightly  negative  (nose  down). 
As  the  LEV  convects  past  the  quarter,  the  pressure  distributions  shifts  so  that  the  upper  surface 
pressure  near  the  leading  is  high  and  the  pressure  downstream  of  the  quarter  chord  is  low, 
causing  the  negative  moment  to  suddenly  decrease  to  a  peak  Cm  =  -0.18  at  a  ~  19°.  After  the 
vortex  has  convected  past  the  airfoil,  for  19°  <  a  <  20°,  Cm  increases  to  -0.12.  As  the 
downstroke  begins,  Cm  gradually  increases  until  the  upper  surface  flow  reattaches  after  which  the 
measured  Cm  is  relatively  weak  and  approximately  constant. 
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Figure  11:  Phase-averaged  CL  (upper  left),  CM  (lower  left),  Cp  contours  (middle)  and  angle 
of  attack  for  M  =  0.2,  a  =  9°  -  11°  cos  (p ,  and  k  =  0.050. 

3.1.2  Dynamic  Response  Considerations 

Temporal  resolution  of  the  dynamic  stall  process  measured  via  the  remote  pressure  sensing 
scheme  outlined  in  Section  2.7  is  one  of  the  ultimate  concerns  which  needs  to  be  addressed 
before  further  analysis  of  the  data.  Assuming  that  it  is  sufficient  to  resolve  content  at  frequencies 
5-6  multiples  of  the  first  harmonic  of  the  pitch  oscillations  (20  Hz),  a  sufficient  frequency 
response  should  be  flat  out  to  approximately  120  Hz.  Furthermore,  the  pneumatic  distortion 
associated  with  the  remote  pressure  sensing  scheme  is  a  strong  function  of  the  mean  pressure 
(Bergh  and  Tijdeman,  1965).  Thus,  the  maximum  distortion  is  expected  for  low  freestream  Mach 
numbers  (generally  high  tunnel  and  model  surface  pressures).  Figure  12  displays  two  oscillation 
periods  of  raw  pressure  traces  recorded  for  the  first  ten  taps  on  the  upper  surface  (including  the 
leading  edge  tap)  for  the  case  of  M=  0.2,  a  =  9°  -  14°  cos  (p,  and  k  =  0.135  (f=  19.5  Hz).  The 
leading  tap  is  designated  Tap  #1  and  tap  number  increases  marching  along  the  upper  surface  in 
the  downstream  direction.  An  offset  equal  to  the  product  of  0.5  psi  and  the  tap  number  is  added 
to  each  pressure  trace  to  spread  out  the  traces  and  improve  visibility.  Traces  are  compared 
without  (left)  and  with  (right)  pneumatic  compensation.  The  difference  between  traces  with  and 
without  pneumatic  compensation  is  subtle.  High  frequency  fluctuations  are  smoothed  over  as  if 
they  were  low-pass  filtered,  and  the  peak  negative  pressures  observed  just  before  the  onset  of 
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stall  (t/T  =  0.5  and  1.5)  are  also  slightly  attenuated.  Figure  13  compares  the  corresponding 
Fourier  series  coefficients  for  the  tap  indicating  the  most  power  of  high-frequency  pressure 
fluctuation  (tap  #2,  see  Appendix  2).  The  upper  plot  in  this  figure  shows  the  Fourier  coefficients 
scaled  by  a  factor  of  five  in  order  to  include  the  compensation  transfer  function  (amplitude  ratio) 
on  the  same  plot.  Evidently  95%  of  the  power  in  this  signal  is  resolved  within  0  </<  200  Hz 
through  the  first  ten  harmonics  of  the  oscillation  frequency.  The  amplitude  ratio  function 
indicates  negative  attenuation  (resonance)  with  a  peak  amplitude  ratio  of  4  at  approximately 
300Hz.  The  amplitude  ratio  is  1.84  at  200Hz  and  the  phase  is  -19°.  This  means  that  high 
frequency  content  is  physically  amplified  in  this  tubing-transducer  configuration.  Whereas  most 
remote  pressure  sensing  schemes  involve  positive  attenuation  and  a  potential  loss  of  high 
frequency  information,  the  frequency  response  in  this  case  is  beneficial  because  high-frequency 
content  is  preserved.  Also,  the  power  in  high  frequencies  is  not  overwhelming  in  the  sense  that  it 
would  distort  lower  frequency  content.  Thus,  the  pressure  traces  in  Figure  12  appear  to  be 
filtered  when  compensation  is  applied.  This  analysis  suggests  that  the  dynamic  response 
characteristics  of  the  tubing-transducer  system  are  adequate  to  resolve  the  unsteady  dynamics  in 
this  flow  scenario. 
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Figure  12:  Comparing  two  periods  of  raw  pressure  traces  for  M=  0.2,  a  ~  9°  -  14°  cos  (p , 

and  k  =  0.135  (f=  19.5  Hz). 
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Figure  13:  Spectra  of  pressure  traces  (tap  #2),  plotted  along  with  pneumatic  compensation 
transfer  function.  M=  0.2,  a  ~  9°  -  14°  cos  (p ,  and  k  =  0.135  (f=  19.5  Hz). 

3.1.3  Validation  with  Historical  Data 

Figure  14  through  Figure  16  are  comparisons  of  OSU  data  and  historical  data  which  has  been 
digitally  extracted  from  graphical  images  in  Lorber  and  Carta  (1987),  also  cited  in  some  figures 
as  AFOSR  TR  87-1202.  Comparisons  are  made  for  M  =  0.2,  a  ~  10°  -  10°  cos  (p,  and  reduced 
frequencies  (k)  equal  to  0.025,  0.050,  and  0.100.  The  angle  of  attack  oscillations  for  k  =  0.025 
and  0.050  have  similar  maximum  a  between  the  OSU  and  Lorber  and  Carta’s  data.  However,  for 
the  case  of  k  =  0.100,  the  maximum  a  of  the  OSU  oscillation  is  2°  higher.  Although  the  OSU 
pitch  oscillation  is  periodic,  it  is  difficult  to  control  the  oscillation  amplitude  because  of 
structural  vibration  of  the  pitch  oscillation  assembly.  For  any  of  the  selected  cases,  there  is  a 
discrepancy  in  CL  which  increases  throughout  the  upstroke.  Also,  the  OSU  data  indicates  a 
slightly  negative  CM ■  during  the  upstroke  whereas  the  Lorber  and  Carta  data  indicates  a 
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C -M  ~  0  .  These  two  differences  are  attributed  to  the  spatial  resolution  of  surface  pressure  near 

the  leading  edge  of  the  OSU  airfoil  model  (see  Section  2.9).  The  onset  of  lift  stall  and  moment 
stall  is  similar  between  both  datasets  for  k  =  0.025  and  0.050.  Generally,  the  OSU  data  indicates 
stall  onset  during  the  upstroke  at  an  a  ~  0.5°  higher  than  indicated  by  the  Lorber  and  Carta  data. 
For  the  case  of  k  =  0.100,  both  datasets  indicate  that  the  flow  is  attached  throughout  the  entire 
upstroke  and  stall  occurs  at  the  maximum  a  as  the  airfoil  begins  the  downstroke.  Apart  from  the 
interval  of  LEV  fonnation  and  convection,  the  post-stall  Cl  and  Cm  values  are  of  similar 
magnitude  during  the  downstroke,  but  the  OSU  data  indicates  lower  moment  magnitudes  in  all 
three  cases  prior  to  reattachment.  During  this  portion  of  the  pitch-oscillation  period,  the 
contributions  to  the  Cl  and  Cm  by  the  pressure  distribution  about  the  leading  edge  are 
diminished.  Thus,  the  error  due  to  the  leading  edge  spatial  resolution  is  less  significant,  and  the 
agreement  between  the  two  datasets  is  improved.  Reattachment  trends  are  also  similar  between 
the  two  datasets.  As  reduced  frequency  increases,  reattachment  generally  occurs  later  during  the 
oscillation  period.  Reattachment  seems  to  occur  slightly  earlier  for  the  OSU  model.  The 
reattachment  angle  of  attack  is  2°,  1°,  and  1°  higher  in  the  OSU  data  for  k  =  0.025,  0.050,  and 
0.100,  respectively.  The  peak  negative  moment  spike  is  of  particular  concern,  and  there  are 
notable  discrepancies  between  the  two  datasets.  In  all  cases,  the  OSU  measurements  indicate 
lower  moment  spike  magnitudes.  The  magnitude  of  the  peak  moment  of  the  OSU  data  is  41%, 
50%,  and  32%  lower  than  the  Lorber  and  Carta  measurement  for  k  =  0.025,  0.050,  and  0.100, 
respectively. 


Figure  14:  Comparing  OSU  data  to  historical  data:  lift  (left)  and  moment  (right)  coefficient 
orbits  for  M=  0.2,  a  ~  10°  -  10°  cos  (p ,  and  k  =  0.025. 
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M  =  0.2  k  =  0.050 


M  =  0.2  k  =  0.050 


Figure  15:  Comparing  OSU  data  to  historical  data:  lift  (left)  and  moment  (right)  coefficient 
orbits  for  M=  0.2,  a  ~  10°  -  10°  cos  (p ,  and  k  =  0.050. 


Figure  16:  Comparing  OSU  data  to  historical  data:  lift  (left)  and  moment  (right)  coefficient 
orbits  for  M  =  0.2,  a  ~  10°  -  10°  cos  (p ,  and  k  =  0.100. 

Figure  17  compares  phase-averaged  pressure  traces  between  both  datasets  for  the  case  of  M  = 
0.2,  a  ~  10°  -  10°  cos  (p,  and  k  =  0.050.  The  relative  validity  of  this  comparison  is  optimal 
because  the  pitch  oscillation  amplitude  and  maximum  angle  are  a  closer  match  than  for  k  = 
0.100.  Traces  are  selected  for  upper  surface  pressure  tap  locations  ( x/c  ~  0.025,  0.220,  0.460) 
which  are  co-located  to  within  0.5%  of  chord  between  the  two  airfoil  models.  Traces  are  offset 
for  visibility.  The  peak  suction  registered  by  each  of  the  pressure  taps  during  the  upstroke  (0  < 
t/T  <  0.5)  is  comparable  between  the  two  datasets.  Flowever,  the  peaks  occurring  in  the  OSU 
measurements  exhibit  a  progressive  phase  lag  relative  to  the  Lorber  and  Carta  measurements. 
The  tap  at  x/c  =  0.025  shows  the  collapse  of  the  suction  peak  near  the  leading  edge  occurring  at 
approximately  the  same  phase  (t/T  =  0.37),  and  the  slope  of  the  pressure  rise  with  time  is  only 
slightly  steeper  for  the  Lorber  and  Carta  data.  The  pressure  rise  signal  convects  past  the 
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subsequent  pressure  taps  of  Lorber  and  Carta  more  quickly  than  observed  in  the  OSU  taps.  This 
is  consistent  with  the  delayed  onset  of  stall  observed  for  the  OSU  data  in  Figure  15.  Both  models 
show  evidence  of  stall  for  0.45  <  t/T  <  0.7,  although  the  Lorber  and  Carta  pressures  are 
consistently  higher  during  this  interval.  This  behavior  explains  why  Lorber  and  Carta’s  moment 
coefficient  is  slightly  more  negative  during  the  beginning  of  the  downstroke.  Earlier 
reattachment  in  the  OSU  data  is  evident  beginning  at  t/T  =  0.62  at  x/c  =  0.025,  and  the 
reattachment  starts  near  the  leading  edge  and  progresses  downstream.  The  difference  between 
the  two  models  in  attached  flow  pressures  is  negligible. 

While  an  unsteady  interference  evaluation  has  not  been  perfonned  for  testing  in  the  6”  x  22” 
transonic  wind  tunnel,  it  is  speculated  that  the  reason  why  suction  levels  on  the  OSU  model  are 
persistently  higher  than  on  Lorber  and  Carta’s  model  is  associated  with  the  aspect  ratio.  The 
aspect  ratio  is  1  and  5.5  for  the  OSU  and  Lorber  and  Carta’s  models,  respectively.  The  low 
aspect  ratio  can  corrupt  the  two-dimensionality  of  the  flow  as  the  airfoil  stalls,  and  can  promote 
artificially  high  velocity  in  the  midspan  region  as  endwall  boundary  layers  effectively  constrict 
the  midspan  flow.  In  the  context  of  the  dynamic  stall  experiment,  this  impact  of  low  aspect  ratio 
may  be  to  delay  the  onset  of  dynamic  stall,  promote  early  dynamic  reattachment,  and  weaken  the 
LEV  vortex,  thereby  reducing  the  strength  of  the  pitching  moment  spike.  A  suitable  interference 
evaluation  might  consist  of  unsteady  testing  using  congruent  airfoil  models  of  increasing  aspect 
ratio  to  detennine  the  aspect  ratio  above  which  the  time -resolved  Cp  is  constant.  For  lack  of 
interference  information,  it  is  sunnised  that  the  measurement  methodology  is  adequate  to  resolve 
the  unsteady  fluid  dynamics,  while  the  scale  of  the  model  may  be  reduced  to  improve  the  validity 
of  future  experiments. 
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Phase-Averaged  Pressure  Traces 


Figure  17:  Phase-averaged  pressure  traces  from  OSU  and  historical  data.  Traces  for  three 
chordwise  positions,  offset  for  visibility.  Test  condition: 

M  =  0.2,  a  ~  10°  -  10°  cos  (p,  and  k  =  0.050. 

3.1.4  Compressibility  Effects 

The  results  acquired  for  the  oscillating  clean  airfoil  in  a  steady  airstream  are  presented  in  two 
ways.  Lift  and  moment  coefficient  orbits  are  presented  in  Figure  18  through  Figure  20  for  fixed 
freestream  Mach  number  to  illustrate  the  impact  of  reduced  frequency.  Results  are  then  presented 
for  fixed  reduced  frequency  (Figure  21  through  Figure  23)  to  understand  how  the  dynamic  stall 
process  is  impacted  by  increasing  the  freestream  Mach  number  and  elevating  the  role  of 
compressibility  effects. 

Dynamic  stall  experiments  were  performed  for  three  values  of  freestream  Mach  number, 
namely  0.2,  0.4,  and  0.6.  In  all  tests  the  target  pitch  oscillation  profile  is  a  ~  9°  -  1 1°  cos  cp,  with 
the  target  amplitude  dictated  by  the  geometry  of  the  cam  in  the  pitch  oscillation  assembly  (see 
Section  2.2).  As  mentioned  previously,  the  maximum  angle  of  attack  is  set  to  20°  with  the  pitch 
oscillation  unpowered.  Attempts  to  minimize  the  vibration  of  the  pitch  oscillation  assembly 
during  operation  were  productive.  However,  vibrations  persist  which  cause  the  measured  angle 
of  attack  amplitude  to  increase  with  frequency.  The  difficulty  in  controlling  the  pitch  amplitude 
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can  make  it  challenging  to  parse  the  effects  of  reduced  frequency.  For  instance,  as  the  pitching 
amplitude  increases,  the  angular  acceleration  and  effective  frequency  at  the  onset  of  stall  or 
reattachment  may  be  slightly  higher.  Generally,  test  dynamic  stall  experiments  were  performed 
in  increments  of  k  =  0.025  up  to  the  limit  of  the  pitch  oscillation  assembly. 

For  M=  0.2  (Figure  18),  data  points  are  presented  for  five  k  values  up  to  a  maximum  of  k  = 
0.125.  As  reduced  frequency  increases,  the  onsets  of  both  dynamic  stall  and  dynamic 
reattachment  are  delayed.  For  k  >  0.075,  attached  flow  is  sustained  up  to  the  maximum  angle  of 
attack  during  the  upstroke,  and  dynamic  stall  occurs  as  the  airfoil  begins  the  downstroke.  As 
expected,  hysteresis  effects  become  more  exaggerated  with  increasing  frequency.  The  post-stall 
Cl  decreases  indicating  a  stronger  stall  and  larger  fluctuations  in  lift  force  throughout  a  pitching 
cycle.  The  peak  negative  pitching  moment  increases  with  reduced  frequency.  For  the  attached 
flow  interval  of  the  oscillation  period,  the  sign  of  the  measured  moment  flips  between  upstroke 
(nose  down)  and  downstroke  and  the  magnitude  increases  with  frequency  during  the  downstroke. 

Figure  19  presents  data  for  three  reduced  frequencies  up  to  a  maximum  of  k  =  0.070  for  M  = 
0.4.  The  onsets  of  stall  and  reattachment  are  similarly  delayed.  Stall  occurs  before  the  airfoil 
reaches  peak  angle  of  attack  in  all  cases  presented.  The  portion  of  either  the  lift  or  moment 
curves  during  the  onset  of  stall  is  referred  to  as  the  stall  characteristic.  For  M  =  0.4,  the  stall 
characteristic  in  either  Cl  or  Cm  is  of  similar  slope  for  all  reduced  frequencies.  As  the  timescale 
of  the  airfoil  motion  is  reduced  relative  to  the  fluid  dynamic  timescales,  the  slope  of  the  moment 
stall  characteristic  decreases  in  magnitude.  Moreover,  the  range  of  peak  negative  moment  values 
is  only  0.035  compared  to  0.07  for  the  same  range  of  k  at  M=  0.2.  As  with  M=  0.2,  the  disparity 
in  CM  values  between  attached  flow  portions  of  the  upstroke  and  downstroke  increases  with 
reduced  frequency.  However,  unlike  the  M  =  0.2  orbits,  there  is  a  notable  difference  in  the  in  Cl 
values  between  attached  flow  portions  of  the  upstroke  and  downstroke  for  M=  0.4.  Generally,  Cl 
is  0.2  higher  after  dynamic  reattachment  than  during  the  upstroke  (a  <  5°),  and  the  dissimilarity 
increases  slightly  with  reduced  frequency. 

Only  two  conditions  are  presented  in  Figure  20  for  M  =  0.6.  Similar  assertions  can  be  made 
about  the  effect  of  increasing  reduced  frequency  regarding  the  delay  of  dynamic  stall  and 
reattachment.  As  with  M  =  0.4,  the  Cl  and  Cm  stall  characteristics  have  similar  slope  for  both 
reduced  frequencies.  Evidently,  enhanced  compressibility  brings  about  a  pronounced  CL 
hysteresis  effect  during  the  attached  interval  of  the  cycle.  The  Cl  measured  after  the  dynamic 
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reattachment  is  consistently  higher  than  the  values  measured  during  the  upstroke.  The  largest 
difference  in  CL  between  upstroke  and  downstroke  is  approximately  0.3  for  k  =  0.025  (a  ~  4°) 
and  0.4  for  k  =  0.050  (a  ~  0°). 

Figure  21,  Figure  22,  and  Figure  23  are  organized  to  show  Cl  and  Cm  hysteresis  loops  for 
varied  M  and  fixed  k  =  0.025,  0.050,  and  0.075,  respectively.  Orbits  of  k  =  0.025  are  inspected 
first.  Increasing  freestream  Mach  number  causes  stall  earlier  during  the  upstroke.  The  maximum 
Cl  occurs  at  a  =  17.5°,  16.5°,  and  14°  for  M  =  0.2,  0.4,  and  0.6,  respectively.  For  a  higher  than 
the  value  of  maximum  Cl,  the  lift-curve  slope  deviates  from  the  pre-stall  value.  On  the  other 
hand,  the  departure  of  CM  from  pre-stall  values  precedes  lift  stall  depending  on  Mach  number. 

For  M  =  0.2,  the  departure  of  CL  and  CM  coincide  at  approximately  17.5°.  For  M  =  0.4,  the 

departure  of  CMa  occurs  approximately  1°  before  that  of  CL  .  For  M  =  0.6,  CMa  changes 

approximately  1.5°  before  CL  .  The  slope  of  both  C/  and  Cm  stall  characteristics  decrease  in 

magnitude  as  M  increases.  The  steep  characteristic  of  the  M  =  0.2  case  is  indicative  of  leading 
edge  stall  behavior  which  is  typically  experienced  by  an  airfoil  executing  sinusoidal  or  ramping 
motion  through  the  stall  boundary  in  an  incompressible  flow.  In  this  flow  scenario,  wall  bound 
vorticity  coalesces  at  the  leading  edge  into  a  strong  spanwise  vortex  (the  LEV).  The  collapse  of 
the  suction  peak  near  the  leading  edge  coincides  with  the  downstream  convection  of  the  LEV, 
which  leaves  behind  an  open  separation  emanating  from  the  leading  edge.  As  compressibility 
effects  become  more  important,  flow  around  the  leading  edge  becomes  supercritical.  If  the  flow 
expansion  dictated  by  the  geometry  and  pressure  gradient  is  sufficiently  strong  a  nonnal  shock 
will  form  on  the  upper  surface  to  facilitate  the  pressure  rise.  Furthennore,  if  the  flow  speeds 
developed  around  the  leading  edge  are  sufficiently  high  (local  M  >  1.2),  the  shock  wave  / 
boundary  layer  interaction  can  induce  stall.  In  the  scenario  involving  an  oscillating  airfoil,  the 
point  of  incipient  stall  does  not  begin  at  the  leading  edge,  though  it  does  gradually  migrate 
upstream  as  a  continues  to  increase.  The  shock  is  strongest  (largest  pressure  rise  across)  and 
furthest  downstream  when  streamline  curvature  around  the  leading  edge  is  maximum.  When  the 
boundary  layer  thickens,  the  effective  geometry  seen  by  oncoming  flow  is  altered  such  that  the 
streamline  curvature  decreases  as  does  the  shock  strength.  The  shock  also  moves  upstream.  Thus, 
as  the  stall  develops  and  grows  in  size,  the  shock  moves  upstream  and  disappears.  The  stall 
behavior  is  tenned  a  trailing  edge  type  since  it  does  not  emanate  from  the  leading  edge  and  it 
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produces  stall  characteristics  which  are  less  abrupt.  The  stall  process  transformation  from 
leading-edge  to  trading-edge  type  can  inhibit  the  formation  and  strength  of  the  LEV  and  result  in 
lower  peak  negative  moment  magnitude  as  evident  in  Figure  21.  The  post-stall  Cl  and  Cm  values 
are  remarkably  similar.  Post-stall  Cl  and  Cm  values  collapse  for  all  three  Mach  numbers  to 
approximately  CL  =  0.8  during  the  downstroke  (20°  >  a  >  10°)  and  CM  =  -0.12  at  a  =  20°.  The 
reattachment  behavior  indicated  in  Cm  are  more  similar  from  M  =  0.4  and  0.6,  with  pitching 
moment  returning  to  attached  flow  values  at  a  ~  10°.  Meanwhile  Cm  returns  to  attached  flow 
values  closer  to  a  ~  12°  for  M  =  0.2.  Perhaps  more  striking  is  that  dynamic  reattachment  is 
evident  in  Cl  at  a  ~  9°  for  all  three  Mach  numbers.  The  effect  of  compressibility  on  the 
hysteresis  behavior  of  the  attached  flow  Cl  is  also  very  clear  even  at  the  lowest  reduced 
frequency. 

Comparisons  are  also  made  in  Figure  22  across  all  three  Mach  numbers  for  increased  k  = 
0.050.  As  with  k  =  0.025,  the  stall  characteristics  exhibit  the  same  morphology  with  increasing 
Mach  number.  For  M=  0.2,  stall  begins  just  before  the  maximum  angle  of  attack  at  a  =  19.5°.  As 
M  increases,  stall  occurs  earlier  in  phase,  and  the  stall  characteristic  flattens  out  to  resemble  the 
trailing  edge  stall  behavior.  The  lift  and  moment  stall  angle  is  again  considered  to  be  the  angle 
during  the  upstroke  at  which  the  slope  of  the  respective  curves  changes  from  the  pre-stall  to  post¬ 
stall  value.  Moment  stall  precedes  lift  stall  by  0.5°  for  M  =  0.2,  1°  for  M=  0.4,  and  2°  for  M  = 
0.6  which  is  a  trend  similar  to  that  of  k  =  0.025.  Post-stall  Cl  during  the  downstroke  collapse  to 
approximately  0.65,  and  peak  negative  moments  are  all  approximately  Cm  =  -0.17  for  k  =  0.050. 
The  reattachment  behavior  is  the  same  as  for  k  =  0.025  as  well.  The  Cl  data  indicates  attached 
flow  Cl  beginning  at  a  —  8°  for  M=  0.2,  and  perhaps  slightly  later  (lower  a)  forM=  0.4  and  0.6. 
The  dynamic  reattachment  of  the  higher  Mach  number  cases  is  more  difficult  to  pinpoint  in  the 
Cl  behavior  because  of  the  attached  flow  hysteresis.  This  figure  most  clearly  shows  how  the  Cl 
disparity  between  upstroke  and  downstroke  during  the  attached  flow  interval  is  a  strong  function 
of  freestream  Mach  number. 

Figure  23  shows  the  hysteresis  loops  with  k  =  0.075  for  M  =  0.2  and  0.4  only  since  the 
maximum  attainable  reduced  frequency  is  k  =  0.050  at  M  =  0.6.  This  is  the  highest  reduced 
frequency  (k  =  0.075)  for  which  valid  comparisons  can  be  made  with  respect  to  the  effect  of 
freestream  Mach  number.  This  figure  continues  to  validate  observed  trends  described  for  k  = 
0.025  and  0.050.  All  generic  assertions  regarding  the  delayed  onset  of  stall,  precedence  of 
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moment  stall,  modified  stall  characteristic,  post-stall  Cl  and  Cm,  and  dynamic  reattachment  are 
consistent  at  this  reduced  frequency. 

Figure  24  through  Figure  29  present  phase-averaged  contours  of  upper  surface  pressure 
coefficient  for  each  Mach  number  (0.2,  0.4,  and  0.6)  at  k  =  0.025  and  0.050  to  illustrate  how  the 
pressure  development  is  modified  with  enhanced  compressibility.  These  figures  are  fonnatted  in 
the  same  way  as  Figure  1 1  with  chordwise  position  along  the  abscissa  and  phase  angle  along  the 
ordinate.  The  contour  intervals  and  colorbars  are  unifonn  for  ease  of  comparison.  Also,  when 
applicable,  a  single  dark  black  contour  line  is  plotted  demarcating  the  critical  pressure  coefficient 
for  a  given  freestream  Mach  number  (Figure  26  -  Figure  29).  The  critical  pressure  coefficient 
corresponds  with  a  local  flow  speed  equal  to  the  speed  of  sound  (M  =  1).  The  critical  pressure 
coefficients  for  M  =  0.2,  0.4,  and  0.6  are  CPiCr  =  -16.3,  -3.7,  and  -1.4,  respectively.  When  the 
local  pressure  is  lower  than  the  critical  pressure,  the  local  Mach  number  over  a  chordwise  extent 
bounded  by  this  contour  is  supercritical. 

Figure  24  and  Figure  25  show  Cp  contours  for  M=  0.2,  k  =  0.025  and  0.050.  In  both  of  these 
cases,  suction  develops  in  the  leading  edge  region  attaining  a  peak  Cp  ~  -7  at  xtc  =  0  while  16°  < 
a  <  17°  during  the  upstroke  just  prior  to  LEV  formation.  As  a  continues  to  increase,  the  peak 
suction  starts  to  diminish,  the  LEV  fonns,  detaches  from  the  leading  edge  and  convects 
downstream  leaving  a  low  pressure  trace  between  140°  <  cp  <  165°.  The  trace  starts  to  fade  as  the 
vortex  convects  past  xtc  =  0.6.  The  LEV  passage  is  immediately  succeeded  by  a  plateau  of 
uniform  Cp  ~  -1  until  the  reattachment  process  begins.  For  k  =  0.025  (Figure  24)  pressure  starts 
to  decrease  first  at  the  leading  edge  at  (p  =  250°  (a  =  14.5°)  and  progresses  toward  the  trailing 
edge  until  (p  =  255°  (a  =  11°)  when  the  flow  is  completely  attached.  The  reattachment  process  is 
similar,  although  delayed  for  k=  0.050. 

Figure  26  and  Figure  27  show  Cp  contours  for  M=  0.4,  k  =  0.025  and  0.050.  Both  cases  have 
a  brief  instance  of  supercritical  flow.  In  both  cases,  the  peak  Cp  ~  -4  at  xtc  =  0.023  while  14°  <  a 
<  16°.  The  depressed  Cp  signal  of  the  LEV  is  evident,  although  like  the  peak  suction,  it  is 
weakened  relative  to  the  M=  0.2  cases. 

Figure  28  and  Figure  29  show  Cp  contours  for  M  =  0.6,  k  =  0.025  and  0.050.  Both  reduced 
frequencies  for  this  Mach  number  feature  significant  intervals  of  supercritical  flow  extending 
over  the  leading  25%  of  chord.  For  k  =  0.025,  the  peak  Cp  ~  -2.5,  and  for  k  =  0.050,  the  peak  Cp 
~  -2.  In  both  cases,  before  the  loss  of  lift  which  occurs  during  the  upstroke,  peak  suction  occurs 
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at  x/c  =  0.074.  Actually,  the  attached  flow  pressure  distribution  near  the  leading  edge  is 
approximated  by  a  plateau  with  Cp  ~  -2  increasing  slightly  from  x/c  =  0.023  to  a  maximum  at  x/c 
=  0.074.  Figure  29  shows  supplemental  upper  surface  Cp  distributions  contrasting  two 
representative  phases  of  supercritical  flow  for  k  =  0.050.  One  of  the  phases  occurs  during  the 
upstroke  (a  =  13°,  (p  =  1 13°)  and  the  other  occurs  during  the  downstroke  (a  =  5.5°,  (p  =  281°).  It 
is  sunnised  that  during  both  the  upstroke  and  downstroke,  the  supercritical  flow  is  tenninated  by 
a  normal  shock  situated  on  the  airfoil  surface  between  0.07  <  x/c  <  0.20  where  the  adverse 
pressure  gradient  is  strong.  During  the  upstroke,  the  shock  gradient  appears  to  be  less  than  that 
which  occurs  during  the  downstroke  and  this  is  evident  by  comparing  the  chordwise  extent  of  the 
supercritical  region  in  the  Cp  contour  plots  for  either  k  =  0.025  or  0.050.  In  the  case  of  both 
reduced  frequencies  suction  diminishes  along  with  the  size  of  the  supercritical  region  beginning 
at  (p  ~  130°  between  15°<a<17°.  However,  this  process  is  much  more  gradual  than  the  lower 
Mach  number  cases  and  there  is  no  trace  of  a  LEV  at  any  phase.  Figure  28  and  Figure  29  also 
clearly  show  the  attached  flow  hysteresis.  For  instance,  in  either  case  ( k  =  0.025  or  0.050)  a  =  8° 
for  cp  =  90°  and  (p  =  270°.  At  <p  =  90°  the  peak  Cp  is  between  -1.6  <  Cp  <  -1.1,  and  at  cp  =  90°  the 
peak  Cp  is  between  -2.3  <  cp  <  -2. 

3.1.5  A  Note  Regarding  Attached  Flow  Hysteresis 

The  attached  flow  hysteresis  is  considered  regarding  whether  the  attached  flow  hysteresis  is  a 
result  of  artificial  post-processing  error  or  some  inadequacy  of  the  experimental  methodology. 
For  instance,  is  it  conceivable  that  the  entire  phase-average  Cp  distribution  is  shifted  in  phase 
because  of  a  measurement  phase  lag  which  is  a  function  of  freestream  Mach  number  and/or  static 
pressure.  The  following  explanation  is  offered  to  support  the  validity  of  the  attached  flow 
hysteresis  phenomenon.  The  low  Mach  number  case  is  consistent  with  patterns  in  historical  data, 
and  this  would  be  the  case  for  which  maximum  signal  distortion  would  be  expected.  The  lift 
curve  collapses  for  the  upstroke  and  downstroke  during  the  attached  flow  interval.  Furthennore, 
the  lift  is  a  strong  function  of  the  integrated  pressure  distribution,  and  a  weak  function  of  angle  of 
attack,  especially  where  a  is  as  small  as  it  is  during  this  interval  (a  <  10°).  Across  all  three  Mach 
numbers  tested,  the  difference  in  attached  flow  lift  as  a  function  of  a  during  the  upstroke  does 
not  exceed  0.1  for  any  k  value.  Comparing  with  the  dissimilarities  documented  above  seems  to 
substantiate  the  claim  of  attached  flow  hysteresis  that  is  a  strong  function  of  compressibility. 
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Figure  18:  M=  0.2,  a  ~  10°  -  10°  cos  <p. 
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Figure  19:  M  =  0.4,  a  ~  10°  -  10°  cos  (p. 
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Figure  20:  M=  0.6,  a  ~  10°  -  10°  cos  <p. 
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Figure  21:  a  ~  10°  -  10°  cos  <p,  and  k  =  0.025. 
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Figure  22:  a  ~  10°  -  10°  cos  <p,  and  k  =  0.050. 
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Figure  23:  a  ~  10°  -  10°  cos  <p,  and  k  =  0.075. 
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Figure  24:  M  =  0.2,  a  =  9°  -  11°  cos  <p ,  and  k  =  0.025. 
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Figure  25:  M=  0.2,  a  =  9°  -  11°  cos  (p ,  and  k  =  0.050. 
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Figure  26:  M=  0.4,  a  =  9°  -  11°  cos  (p,  and  k  =  0.025. 


M  =  0.4  a  =9°+/- 13°  k  =  0.050 


x/c  a 

Figure  27:  M=  0.4,  a  =  9°  -  13°  cos  (p,  and  k  =  0.050. 
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M  =  0.6  a  =8°+/- 12°  k  =  0.025 


Figure  28:  M  =  0.6,  a  =  8°  -  12°  cos  <p ,  and  k  =  0.025. 


M  =  0.6  a  =8°+/- 15°  k  =  0.050 


Figure  29:  M=  0.6,  a  =  8°  -  15°  cos  (p ,  and  k  =  0.050. 
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3.2  Mach  Oscillations  with  fixed  Angle  of  Attack 

The  conventional  model  for  understanding  compressible  dynamic  stall  involves  studying  the 
coupled  fluid  and  structural  dynamics  of  an  airfoil  executing  a  prescribed  sinusoidal  oscillation 
of  angle  of  attack  in  a  steady  airstream.  Elucidating  the  impact  of  time-varying  compressibility 
on  compressible  dynamic  stall  ultimately  requires  modeling  a  sinusoidal  airfoil  motion  which  is 
synchronized  with  a  sinusoidal  freestream  Mach  number  oscillation.  Understanding  the 
complexity  of  coupled  airfoil  and  freestream  unsteadiness  is  aided  by  isolating  the  effects  of  each 
aspect.  Therefore,  before  presenting  results  with  synchronized  airfoil  pitching  motion  and 
freestream  Mach  number  oscillation,  this  section  focuses  on  a  fixed-pitch  airfoil  subjected  to  the 
time-varying  airstream  condition.  The  data  will  show  that  an  important  angle  of  attack  is  the 
static  stall  angle  for  the  time-mean  Mach  number.  For  a  steady  flow  Mach  number  of  0.4,  the 
static  stall  angle  of  attack  is  measured  to  be  14°  <  ass  <  15°  using  the  clean  SSC-A09  model  in 
the  6”  x  22”  transonic  wind  tunnel.  Oscillating  airstream  data  were  acquired  for  fixed  airfoil 
incidence  angles  ranging  from  pre-static-stall  a  to  post-static-stall  a  to  identify  the  effects  of  the 
time-varying  freestream  with  respect  to  the  static  stall  angle  for  the  mean  Mach  number  of  0.4. 
The  set  of  angles  tested  includes  a  =  10°,  13°,  14°,  15°,  16°,  and  18°.  Given  the  maximum  Mach 
number  oscillation  frequency  of  17  Hz,  unsteady  testing  was  perfonned  for  reduced  frequencies 
of  k  =  0.0125,  0.0250,  and  0.0500.  Phase-averaged  lift  and  moment  coefficient  measurements  are 
presented  along  with  selected  Cp  contour  plots  in  order  to  illustrate  the  impact  of  freestream 
modulation  on  airfoil  perfonnance. 

Figure  30  and  Figure  31  combine  angle  of  attack,  oscillating  Mach  number,  reduced 
frequency,  and  phase-averaged  performance  metric  information  to  present  representative  cases 
with  respect  to  the  static  stall  angle  for  M=  0.4.  In  both  figures,  a  is  plotted  along  the  x-axis  and 
Mis  plotted  along  the  y-axis.  Cl  is  plotted  along  the  z-axis  in  Figure  30,  and  Cm  is  plotted  along 
the  z-axis  in  Figure  31.  The  four  incidence  angles  represented  are  a  =  10°,  13°,  15°,  and  18°. 
Curves  are  color-coded  for  reduced  frequency  with  red  curves  corresponding  to  k  =  0.025  and 
black  curves  corresponding  to  k  =  0.050.  Arrows  are  added  to  the  curves  to  show  the  direction  of 
increasing  phase  angle.  The  flow  is  decelerating  from  maximum  to  minimum  M  for  0°  <  cp  < 
180°  and  accelerating  to  maximum  M  for  180°  <  cp  <  360°.  The  Mach  number  waveform  (M  = 
0.4  +  0.07  cos  (p)  is  presented  as  a  function  of  phase  in  subsequent  figures  and  will  be  discussed 
in  more  detail  later. 


J.  W.  Gregory,  W91 1NF1 1 10503  Final  Report,  p.  51 


Both  figures  show  that  the  hysteretic  behavior,  which  depends  on  freestream  acceleration,  is 
transformed  as  a  is  swept  through  the  static  stall  angle  for  M  =  0.4.  Generally,  both  Cl  and  Cm 
tend  to  oscillate  about  the  steady  flow  values  and  both  are  measured  to  be  greater  while  the 
freestream  is  decelerating  than  while  accelerating.  Especially  for  a  <  ass,  the  hysteresis  is 
exaggerated  by  increasing  the  acceleration  magnitudes  which  are  proportional  to  reduced 
frequency.  Also,  unsteadiness  becomes  strong  enough  for  a  >  ass  to  be  evident  in  the  phase- 
averaged  orbits.  Looking  specifically  at  the  C±  behavior  in  Figure  30,  hysteresis  is  biased  toward 
the  minimum  Mach  number  where  the  Cl  disparity  is  greatest.  For  a  =  15°,  the  performance  is 
straddling  the  stall  boundary.  Thus  the  unsteadiness  is  elevated,  and  the  performance  toggles 
between  pre-stall  values  at  low  M  and  post-stall  values  at  high  M.  This  will  be  discussed  in  more 
detail  later.  For  a  =  18°,  the  performance  oscillates  about  the  post-static-stall  value,  and  the 
hysteresis  is  still  comparable  to  pre-stall  (a  =  13°)  though  not  as  clearly  a  strong  function  of 
reduced  frequency.  The  CM  behavior  portrayed  in  Figure  3 1  shows  that  the  CM  hysteresis  is 
greatest  for  a  =  15°  where  the  performance  straddles  the  stall  boundary,  and  it  is  again  biased 
toward  the  minimum  Mach  number.  Pre-stall  hysteresis  is  evident,  but  relatively  weak,  and  post¬ 
stall  hysteresis  appears  to  be  negligible. 


Figure  30:  Cl  vs.  a  vs.  M  for  representative  cases  of  fixed  alpha  and  oscillating  airstream. 
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Figure  31:  Cm  vs.  a  vs.  M  for  representative  cases  of  fixed  alpha  and  oscillating  airstream. 

Figure  32  and  Figure  33  present  representative  phase-average  contours  of  the  upper  surface 
Cp  distribution  in  a  format  similar  to  plots  presented  in  previous  sections.  Chordwise  position  is 
plotted  along  the  abscissa  and  phase  angle  along  the  ordinate.  The  phase-averaged  Mach  number 
is  plotted  to  the  right  with  Mach  number  along  the  abscissa  and  phase  angle  along  the  ordinate. 
As  mentioned  earlier,  the  phase  angle  is  ordered  such  that  the  freestream  flow  is  decelerating  for 
0°  <  (p  <  180°  and  accelerating  for  180°  <  cp  <  360°.  The  contour  levels  and  colorbars  are 
consistent  between  these  two  figures  for  ease  of  comparison.  Figure  32  presents  results  for  a 
selected  pre-stall  angle  of  attack,  a  =  14°,  with  k  =  0.025.  Figure  33  presents  a  data  point  for 
which  the  performance  straddles  the  static  stall  boundary  (a  =  15°)  and  k  =  0.025.  Although  the 
surface  pressures  do  fluctuate  slightly  for  post-stall  a,  the  upper  surface  Cp  distribution  plateaus 
across  the  chord  at  approximately  Cp  ~  -1  indicating  an  open  separation. 


J.  W.  Gregory,  W91 1NF1 1 10503  Final  Report,  p.  53 


a  =  14°  k  =  0.0261  f=  7.1Hz 


Figure  32:  Contours  of  upper  surface  Cp  for  a  =  14 °,M=  0.4  ±  0.07,  and  k  =  0.025. 


Figure  33:  Contours  of  upper  surface  Cp  for  a  =  15°,  M  =  0.4  ±  0.07,  and  k  =  0.025. 
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The  pattern  displayed  in  Figure  32  for  a  =  14°  is  typical  for  pre-stall  a  in  the  sense  that  the 
greatest  suction  peak  is  measured  to  be  approximately  Cp  ~  -5.5  at  xtc  =  0.023  during  the 
deceleration  interval  between  90°  <  cp  <  180°.  The  maximum  magnitude  of  freestream 
deceleration  occurs  at  approximately  cp  ~  90°.  As  the  flow  accelerates,  peak  suction  decreases  as 
do  Cl  and  Cm-  This  angle  of  attack  (a  =  14°)  was  also  selected  because  it  shows  the  signs  of 
unsteadiness  surviving  the  averaging  process  as  Mach  number  approaches  the  maximum  value. 
Specifically,  when  unsteadiness  signals  through  the  averaging  process  it  produces  striations  in 
the  contours  which  stretch  along  chordwise  direction  and  fluctuate  in  phase. 

As  a  increases  to  15°,  the  unsteadiness  is  amplified  producing  the  wavy  contours  in  Figure 
33.  Similar  to  pre-stall  a,  maximum  suction  develops  near  the  leading  edge  during  the 
deceleration  interval.  However,  unlike  the  gradual  loss  of  suction  observed  during  the 
acceleration  interval  for  pre-stall  a,  acceleration  at  a  =  15°  causes  the  suction  peak  to  decrease 
suddenly  at  cp  =  235°,  when  M  ~  0.36  followed  by  fluctuating  pressure  distribution  near  the 
leading  edge. 

The  airfoil  response  to  freestream  modulation  for  the  tested  range  of  a  at  k  =  0.025  is 
condensed  further  in  Figure  34.  Four  plots  are  stacked  so  that  the  phase  correlation  can  be  easily 
interpreted.  The  lowest  plot  is  phase-averaged  Mach  number  versus  phase  angle.  Above  Mach 
number,  are  CM  and  CL.  Finally,  the  uppermost  figure  shows  the  phase-averaged  Cp  measured  by 
the  upper  surface  pressure  tap  located  at  xtc  =  0.074.  The  Cp  curves  show  generally  lower 
pressures  during  deceleration.  As  a  increases  from  10°  to  13°  there  is  a  bulk  negative  shift  in 
suction  and  a  positive  shift  in  lift  which  is  expected  for  pre-stall  alpha.  As  a  increases  from  13° 
to  14°,  suction  and  lift  are  approximately  the  same  for  M  <  0.42,  but  decreased  for  M  >  0.42. 
This  trend  continues  for  a  increased  to  15°  as  the  decrease  in  suction  and  lift  is  even  more 
dramatic  for  M  >  0.4  (0°  <  tp  <  90°  and  225°  <  <p  <  360°).  Furthermore,  the  evidence  of 
unsteadiness  is  pronounced  while  M>  0.4.  For  post-stall  a,  the  1/period  pressure  fluctuations  are 
reduced  and  unsteadiness  persists.  The  Cm  behavior  is  readily  divided  into  the  three  distinct 
categories  of  pre-stall,  stall,  and  post-stall  angles.  For  a  <  ass,  Cm  is  approximately  -0.01  with 
relatively  small  1/period  fluctuation.  For  a  >  ass.  Cm  is  approximately  equal  to  the  post-static  stall 
value.  The  moment  coefficient  is  approximately  -0.09  for  a  =  16°  and  Cm  ~  -0.11  for  a  =  18°. 
As  with  the  pre-stall  a,  the  1/period  fluctuations  are  negligible.  The  most  dramatic  1/period 
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fluctuations  of  Cp,  CL  and  CM  are  measured  for  a  =  15°  when  Mach  number  crosses  0.4  (either 

accelerating  or  decelerating)  as  the  airfoil  performance  straddles  the  stall  boundary. 

Phase-Averaged  Cp  at  x/c  =  0.074  for  k  =  0.025 


Phase-Averaged  CM  for  k  =  0.025 

1  "  1  t  E  1  1  1  "  E  "  1  t  E  1  1  E  t  1  "  F 


a  =  10°  a  =  13°  a  =  14°  a  =  15°  a  =  16°  a  =  18° 


Phase-Averaged  M  for  k  =  0.025 


a  =  10° 

a  =  13° 

a  =  14°  . 

.  a  =  15° 

a  =  16° 

a  =  1 8° 

Figure  34:  Phase-averaged  Cp  traces  ( x/c  =  0.074),  Cl,  Cm,  and  M  vs.  phase  angle  for 
M  =  0.4  ±  0.07  and  k  =  0.025  and  representative  a. 


J.  W.  Gregory,  W91 1NF1 1 10503  Final  Report,  p.  56 


The  fluctuating  performance  metrics  are  quantified  in  Figure  35,  where  the  time-averaged  Cl 
and  Cm  are  plotted  with  the  steady  M  =  0.4  clean  airfoil  results.  Error  bars  are  used  to  indicate 
the  range  of  measured  Cl  and  Cm  values  for  each  reduced  frequency  ( k  =  0.0250  and  0.0500). 
For  each  fixed  k  value,  the  Cl  range  function  is  similar  pre-stall  and  post-stall  a.  The  maximum 
range  occurs  at  fixed  a  within  the  stall  characteristic.  The  range  widens  as  reduced  frequency 
increases,  and  the  upper  bound  which  is  generally  measured  during  deceleration,  is  modified  the 
most.  The  Cm  range  is  smallest  for  pre-stall  a.  When  a  straddles  the  static  stall  boundary  with  a  = 
15°  the  phase-averaged  Cm  range  is  maximum  at  A  Cm  ~  0.05  for  k  =  0.025  and  A  Cm  ~  0.072  for 
k  =  0.050.  This  is  because  the  pressure  distribution  is  toggling  between  conditions  of  attached 
and  separated  flow  near  the  leading  edge  producing  large  variations  in  pitching  moment  at  this 
angle. 


(a) 
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Figure  35:  Comparison  of  range  of  measured  (a)  Cl  and  (b)  Cm  for  fixed  a  and 
M  ~  0.4  ±  0.07  and  varied  reduced  frequency;  steady  M=  0.4  results  also  shown. 
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3.3  Coupled  Pitch  and  Mach  Oscillations 

The  final  stage  of  experiments  performed  with  the  clean  SSC-A09  model  is  that  of  airfoil 
pitching  motion  synchronized  with  freestream  oscillation  to  elucidate  the  impact  of  the  time- 
varying  compressibility  condition  on  compressible  dynamic  stall.  All  of  the  results  presented  for 
an  oscillating  airfoil  in  a  steady  airstream  were  intended  to  simulate  the  dynamic  stall 
experienced  by  an  airfoil  executing  a  sinusoidal  pitch  oscillation  with  angle  of  attack 
approximated  by  a  =  8.5°  -  1 1.5°  cos  tp.  All  of  the  results  presented  for  a  fixed  pitch  airfoil  were 
subjected  to  a  Mach  number  oscillation  approximated  by  M=  0.4  +  0.08  cos  (p.  By  this  definition 
of  the  phase  angle,  tp,  the  maximum  angle  of  attack  and  maximum  Mach  number  are  180°  out  of 
phase.  Figure  36  depicts  the  scenario  for  which  this  phase  angle  convention  is  defined.  A 
rotorcraft  is  considered  to  be  flying  in  straight  forward  flight,  such  that  the  relative  velocity  seen 
by  an  airfoil  section  on  the  advancing  blade  is  augmented  by  the  forward  flight  speed,  and  the 
relative  velocity  seen  by  the  retreating  blade  is  reduced.  The  azimuthal  coordinate  is  selected  so 
that  (p  =  0°  corresponds  to  the  blade’s  position  where  the  relative  Mach  number  is  maximum,  and 
tp  =  180°  corresponds  the  minimum  relative  Mach  number.  A  phase  shift  (A (p)  between  the  pitch 
oscillation  and  Mach  number  oscillation  is  defined  such  that  A (p  =  0°  describes  the  condition 
wherein  the  minimum  angle  of  attack  coincides  with  the  maximum  relative  Mach  number  and 
vice  versa.  Generally,  this  was  the  targeted  phase  relationship  for  the  coupled  oscillation 
experiments.  However,  obtaining  the  desired  phasing  during  testing  proved  challenging.  As 
outlined  in  Section  2.4,  the  servo  motor  driving  the  Mach  number  oscillation  assembly  is 
enslaved  to  follow  the  spindle  assembly  of  the  pitch  oscillation  device.  The  phasing  is  adjusted 
prior  to  testing  with  both  units  unpowered.  The  phase  shift  between  the  test-section  Mach 
number  waveform  and  the  choke  vane  position  is  known  from  empty-tunnel  tests  (see  Gompertz 
et  al.  2011),  and  this  information  was  used  to  inform  the  pre-test  phasing  adjustment.  Still,  the 
target  phasing  was  rarely  achieved,  and  more  development  is  required  to  improve  the  precision 
of  this  aspect  of  synchronization. 

In  addition  to  test  points  acquired  while  targeting  A (p  =  0°,  A (p  was  also  deliberately  set  to  a 
wider  range  of  values  to  understand  the  sensitivity  to  and  impact  of  this  variable.  Figure  37 
illustrates  how  the  pitch  and  Mach  number  oscillations  are  aligned  and  how  the  sign  of  the  phase 
shift  is  defined.  Three  angle  of  attack  profiles  are  shown  (red  curves)  along  with  one  Mach 
number  wavefonn  (black  curve).  The  a  curve  for  A^>  =  0°  depicts  the  nominal  condition.  The 
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phase  shift  is  defined  to  be  positive  when  the  airfoil  oscillation  is  leading  the  relative  Mach 
number  oscillation.  For  instance,  the  maximum  a  leads  the  minimum  Mach  number.  Conversely, 
Atp  <  0  when  the  airfoil  oscillation  lags  the  relative  Mach  number  oscillation.  While  this  figure  is 
useful  when  explaining  how  phasing  is  defined  for  the  results  which  follow,  presenting  phase- 
averaged  Cp  contours  in  this  frame  of  reference  may  not  be  intuitive.  Instead,  for  the  purpose  of 
producing  plots  which  are  of  consistent  formation  and  easy  to  interpret  a  phase  shift  is  applied  to 
align  the  pitch  oscillation  with  a  =  ao  -  a\  cos  <p.  The  same  shift  is  applied  to  the  measured  Mach 
number  oscillation,  so  the  sign  convention  for  A (p  has  the  same  meaning. 


Figure  36:  Diagram  of  rotor  disc,  relative  velocity,  and  adopted  phase  angle  convention. 


Figure  37:  Example  phasing  of  synchronous  sinusoidal  a  (red)  and  M  (black)  oscillations. 
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Coupled  oscillation  testing  was  performed  for  a  prescribed  a  =  8.5°  ±  11.5°  and  M  =  0.4  ± 
0.08.  Three  reduced  frequencies  were  tested  including  k  =  0.0125,  0.0250,  and  0.0500.  Lift  and 
moment  coefficient  orbits  are  compiled  for  each  frequency  in  Figure  38,  Figure  40,  and  Figure 
43.  Corresponding  phase-averaged  upper  surface  Cp  contours  are  presented  for  k  =  0.0125, 
0.0250,  and  0.050  in  Figure  39,  Figure  41,  and  Figure  44,  respectively.  Pressure  coefficient 
contours  are  formatted  consistently  throughout  this  document,  and  all  contour  levels  and 
colorbars  are  consistent  within  this  section.  In  some  cases,  data  are  available  for  numerous  phase 
relationships,  and  representative  cases  are  selected  for  presentation  to  avoid  clutter.  On  the  right- 
hand-side  of  each  of  these  figures  is  a  double  axis  plot  which  shows  the  phase  alignment  between 
the  pitch  oscillation  profile  (dashed  red  line)  and  the  Mach  number  oscillation  for  each  of  the 
selected  cases  with  color-coding  to  match  the  curves  in  Cl  and  Cm  plots. 

Hysteresis  loops  are  presented  for  the  lowest  frequency  ( k  =  0.0125)  in  Figure  38. 
Comparable  dynamic  stall  data  for  the  pitching  airfoil  in  a  steady  M  =  0.4  airstream  is  not 
available  for  this  reduced  frequency.  One  of  the  three  cases  presented  is  a  close  match  with  the 
targeted  phase  relationship  with  the  airfoil  oscillation  lagging  the  freestream  oscillation  by  A cp  = 
-7.5°.  In  another  case,  a  lags  Mby  =  -26.5°.  In  the  third  case,  a  leads  Mby  A  cp  =  +127.3°.  In 
the  first  two  cases  for  which  a  lags  M,  the  freestream  deceleration  occurs  during  the  upstroke, 
and  the  freestream  accelerates  during  the  downstream.  The  dynamic  stall  and  reattachment 
processes  for  these  two  cases  are  very  similar  with  Clmox  ~  1.4  and  stall  onset  at  a  =  16.5°  during 
the  upstroke,  and  reattachment  evident  at  a  =  10°.  Moreover,  the  difference  in  moment  stall 
behavior  is  negligible.  Moment  stall  occurs  at  a  =  16.5°  as  well,  and  the  peak  nose  down 
moments  are  equal  at  Cm  =  -0.14.  The  similarity  between  these  two  cases  despite  the  20° 
difference  in  phase  lag  speaks  to  the  sensitivity  of  the  dynamic  stall  process  to  the  precise  phase 
relationship  between  the  synchronized  pitch  and  freestream  oscillations.  These  cases  show  that 
for  k  =  0.0125,  the  inability  to  match  the  targeted  phase  relationship  does  not  invalidate  the 
results.  The  phasing  represented  by  the  third  case  is  not  merely  academic  in  nature.  Such  a  phase 
relationship  might  provide  insight  into  the  physics  of  a  maneuvering  rotorcraft.  In  the  case  where 
maximum  a  occurs  near  the  peak  relative  Mach  number,  more  lift  may  be  generated  by  the 
advancing  blade  causing  the  rotorcraft  to  bank  left.  In  this  third  case,  the  freestream  is 
accelerating  during  most  of  the  upstroke,  and  decelerating  during  the  downstroke.  Also,  a  peaks 
when  M  is  nearly  maximum.  The  resulting  dynamic  stall  is  in  stark  contrast  with  the  first  two 
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lagging  cases.  The  peak  C/  is  decreased  by  approximately  15%,  the  stall  characteristic  indicates 
a  much  more  gradual  loss  of  lift  beginning  at  a  lower  a  during  the  upstroke.  The  lift-curve 
departs  from  the  attached  flow  trajectory  at  a  =  13°,  flattens  out,  and  lift  starts  to  decrease  at  a  = 
15.5°  with  a  shallow  slope  compared  to  the  first  two  cases.  Post-stall  lift  is  marginally  higher 
than  the  measurements  for  the  first  two  cases.  Moment  stall  occurs  at  a  =  15.5°,  and  the  moment 
stall  characteristic  is  parallel  to  those  of  the  first  two  cases.  In  fact,  all  three  cases  exhibit  equal 
magnitudes  of  maximum  nose  down  moment.  The  Cp  contours  corresponding  to  these  k  =  0.0125 
cases  are  presented  in  Figure  39.  The  first  two  lagging  cases  are  virtually  the  same  with  the 
prominent  depression  signal  of  the  leading  edge  vortex  at  cp  ~  135°.  The  LEV  signature  is  less 
discernable  in  the  third  case  (A^  =  +127.3°).  Also,  the  post  stall  Cp  plateaus  at  a  slightly  lower 
pressure  coefficient  in  accord  with  the  greater  post-stall  Cl  observed  for  this  case. 

There  are  two  important  factors  at  play,  and  it  is  difficult  to  identify  the  role  of  each  at  this 
stage.  One  factor  is  compressibility.  The  fact  that  for  a  >  10°  compressibility  effects  are  lessened 
for  the  first  two  cases  and  greater  for  the  third  cases  can  explain  disparate  lift  generation.  The 
second  factor  is  the  impact  of  freestream  acceleration  on  the  hysteresis  as  in  the  manner  observed 
for  the  fixed-pitch  airfoil  in  the  oscillating  airstream  (Section  3.2).  Although  the  acceleration 
magnitudes  are  a  relative  minimum  for  the  freestream  oscillation  at  this  reduced  frequency,  the 
evidence  of  lift  enhancement  in  the  first  two  cases  is  consistent  with  that  observed  for  a  fixed 
pitch  airfoil  subjected  to  a  decelerating  freestream.  The  post  stall  lift  is  slightly  higher  during  the 
downstroke,  and  this  is  likely  to  be  a  freestream  deceleration  effect.  It  was  shown  in  Section  3.1 
that  for  all  reduced  frequencies,  the  post-stall  lift  coefficient  collapsed  for  the  airfoil  oscillating 
in  pitch  in  a  steady  airstream  with  M  =  0.2,  0.4,  0.6.  Thus,  the  discrepancy  in  post-stall  lift  is  not 
attributed  to  the  instantaneous  freestream  Mach  number.  Furthermore,  results  for  the  fixed-pitch 
airfoil  subjected  to  freestream  oscillation  indicate  that  significant  Cl  hysteresis  is  observed  even 
for  post-stall  conditions. 
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Figure  38:  Cl  and  Cm  orbits  for  a  =  8.5°  ±  11°,  M  =  0.4  ±  0.07,  k  =  0.0125,  and 
representative  phasing  between  synchronous  oscillations. 
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Figure  39:  Cp  contours  for  a  =  8.5°  +  11°,  M  =  0.4  ±  0.07,  k  =  0.0125,  and  representative 

phasing  between  synchronous  oscillations. 

Figure  40  displays  the  CL  and  CM  orbits  for  test  points  with  representative  phasing  and  k  = 
0.0250.  Each  plot  also  includes  comparable  dynamic  stall  data  for  k  =  0.0250  and  the  identical 
airfoil  pitch  oscillation  profile  in  a  steady  M  =  0.4  airstream  (dashed  red  curves).  The  closest 
match  with  the  targeted  phase  relationship  is  the  first  case  in  which  the  airfoil  oscillation  lags  the 
oscillating  Mach  number  by  A cp  =  -28.2°.  During  the  upstroke,  as  a  approaches  stall  onset,  the 
measured  Cl  is  slightly  higher  than  the  steady  M  =  0.4  value.  Again,  in  Section  3.1  it  was  shown 
that  the  lift  curves  collapsed  for  M  =  0.2,  0.4,  and  0.6  during  the  attached  flow  portion  of  the 
upstroke,  so  this  is  likely  a  freestream  deceleration  effect.  The  peak  Cl  is  15%  higher  than  the 
steady  M=  0.4  CLmax,  and  it  occurs  when  a  =  17.25°,  0.75°  higher  than  steady  M=  0.4.  The  lift 
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stall  characteristic  is  steepened.  Post-stall  Cl  and  the  reattachment  during  the  downstroke  are 
virtually  the  same  as  the  steady  M  =  0.4  case.  Moment  stall  is  similarly  delayed  by  A  a  =  0.75° 
relative  to  the  steady  M=  0.4  case,  but  the  rest  of  the  cycle  seems  to  be  unaltered. 

The  second  case  presented  in  this  figure  is  for  a  leading  Mach  number  by  Atp  =  +46.4°.  This 
phase  shift  is  approximately  80°  from  the  first  coupled  case  discussed  above,  but  freestream 
deceleration  is  still  underway  during  the  airfoil  upstroke.  Moreover,  the  maximum  deceleration 
magnitude,  which  occurs  as  Mach  number  crosses  its  mean  value,  coincides  with  the  onset  of 
dynamic  stall  at  a  =  16.5°.  In  the  first  case  (Atp  =  -28.2°),  the  onset  of  stall  occurs  when  the 
Mach  number  is  close  to  the  minimum  and  the  acceleration  parameter  is  approaching  zero  before 
changing  sign.  Therefore,  comparing  these  two  specific  cases  provides  insight  into  the  roles  of 
compressibility  and  acceleration  with  respect  modifying  the  dynamic  stall  process.  As  a 
approaches  stall  onset  during  the  upstroke,  a  slightly  higher  Cl  is  measured  in  both  cases  with 
respect  to  the  steady  M  =  0.4  data.  However,  the  higher  peak  Cl  is  observed  for  the  first  case 
where  the  minimum  Mach  number  coincides  with  the  onset  of  stall.  The  lift-curve  slope  changes 
sign  for  both  cases  at  approximately  a  ~  17°.  The  slope  of  the  lift  stall  characteristic  for  the 
second  case  is  shallower  for  the  first,  and  the  post-stall  Cl  for  19°  <  a  <  20°  at  the  end  of  the 
upstroke  is  higher  by  approximately  0.1.  Otherwise  for  the  remainder  of  the  pitch  oscillation,  the 
difference  in  Cl  between  either  case  and  the  steady  M  =  0.4  case  is  negligible.  As  with  k  = 
0.0125,  all  Cm  orbits  register  equal  peak  negative  moments.  For  k  =  0.0250,  the  peak  negative  Cm 
~  -0.16,  and  for  these  first  two  cases,  the  slope  of  the  moment  stall  characteristic  is 
approximately  the  same  as  the  steady  M  =  0.4  case.  Moment  stall  for  the  second  case  ( Atp  = 
+46.4)  occurs  at  the  same  phase  and  a  =  16.5°  as  for  the  steady  M  =  0.4  case,  while  it  is  delayed 
by  A  a  ~  0.75°  for  the  first  case  (Atp  =  -28.2°). 

Phase-averaged  upper  surface  Cp  contours  are  presented  for  the  k  =  0.025  cases  in  Figure  41. 
The  steady  M  =  0.4  data  is  presented  on  the  far  left  with  coupled  oscillation  cases  in  order  of 
increasing  Atp  positioned  from  left  to  right.  The  first  two  coupled  oscillation  cases  indicate 
enhanced  suction  developing  near  the  leading  edge  and  culminating  with  a  stronger  LEV 
signature  relative  to  the  steady  M=  0.4  LEV.  Suction  is  sustained  until  higher  a  for  Atp  =  -28.2° 
thereby  delaying  moment  stall  and  the  strength  of  the  LEV  appears  to  decay  more  slowly  as  it 
convects  downstream  in  this  case. 
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The  third  case  representing  coupled  oscillations  at  k  =  0.0250  is  for  a  leading  oscillating 
Mach  number  by  Atp  =  +145.4°.  Like  the  third  case  for  k  =  0.125  (A <p  =  +127.3),  the  freestream  is 
accelerating  during  the  upstroke,  decelerating  during  the  downstroke  so  that  a  peaks  when  M  is 
nearly  maximum.  The  hysteresis  is  diminished  for  this  phasing.  The  area  enclosed  by  the  Cl  and 
CM  orbits  reduced.  The  difference  in  CLmax  between  upstroke  and  downstroke  decreases  from 
upwards  of  0.6  for  the  first  two  cases  to  0.3  for  this  third  case.  The  lift  stall  process  initiates  at  a 
~  15°  where  the  lift  curve  departs  from  the  first  two  cases.  The  stall  characteristic  indicates  stall 
which  develops  more  gradually,  and  the  post-stall  Cl  during  the  downstroke  is  higher  than  the 
first  two  cases  and  the  steady  M=  0.4  data.  Moment  stall  occurs  0.75°  before  the  steady  M=  0.4 
case  and  the  slope  of  the  moment  stall  characteristic  is  shallower.  The  first  two  cases  and  the 
steady  M  =  0.4  cases  all  have  a  spike  with  peak  nose  down  Cm  =  -0.17  at  a  ~  18°.  The  spike 
occurs  while  the  LEV  convects  between  the  quarter-chord  and  the  trailing  edge.  Once  the  LEV 
convects  past  the  trailing  edge,  the  moment  coefficient  increases  to  CM  =  -0.15  by  the  phase 
where  a  =  20°  as  the  airfoil  completes  the  upstroke.  Whereas  the  other  cases  all  exhibit  the  Cm 
spike,  the  third  case  (A cp  =  +145.4°)  does  not.  The  moment  coefficient  decreases  to  Cm  =  -0.14 
without  the  higher  frequency  fluctuations  that  accompany  the  LEV  convection  sequence.  The 
reattachment  process  is  similar  for  all  four  cases. 

The  comparison  between  Cp  contours  for  A cp  =  +145.4°  and  the  three  other  data  points  in 
Figure  41  is  very  similar  in  nature  to  that  for  k  =  0.0125  (Figure  39).  Suction  development  is 
limited  near  the  leading  edge,  and  the  LEV  signature  is  indiscernible.  Furthermore  the  Cp  levels 
are  slightly  lower  during  the  stalled  interval  relative  to  the  other  three  cases. 
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Figure  40:  Cl  and  Cm  orbits  for  a  =  8.5°  ±  11.5 °,M  =  0.4  ±  0.07,  k  =  0.0250,  and 
representative  phasing  between  synchronous  oscillations. 
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Figure  41:  Cp  contours  for  a  =  8.5°  ±  11.5°,  M=  0.4  ±  0.07,  k  =  0.0250,  and  representative 

phasing  between  synchronous  oscillations. 

Figure  42  represents  phase-averaged  upper  surface  Cp  data  recorded  by  16  taps  near  the 
leading  edge  for  chordwise  positions  spanning  0.023  <  x/c  <  0.335.  The  figure  consists  of  five 
plots  stacked  for  ease  of  comparison.  Phase  angle  is  plotted  along  the  abscissa.  Phase-averaged 
Mach  number  is  reproduced  in  the  lowest  plot  for  reference  regarding  the  phasing  (A (p)  in  each 
case.  Pressure  coefficient  data  are  presented  for  each  of  the  k  =  0.0250  cases  (one  A  cp  case  in 
each  of  the  upper  three  Cp  plots)  and  for  the  steady  M  =  0.4  case  (lower  Cp  plot).  No  artificial 
offset  is  used  to  add  space  between  curves  of  consecutive  taps  for  clarity  because  the  levels  are 
compared  with  the  phase-averaged  critical  pressure  coefficient  to  ascertain  the  extent  of 
supercritical  flow.  Finally,  it  should  be  noted  that  the  upstream -most  tap  always  registers  the 
lowest  Cp  for  90°  <  (p  <  315°,  and  the  peak  suction  indicated  by  each  tap  decreases  marching  in 
the  downstream  direction.  Attention  is  called  to  the  development  of  suction  with  respect  to  CPjCr 
and  onset  of  stall  occurring  for  45°  <  cp  <  150°.  Two  upper  surface  taps  measure  supercritical 
pressures  for  108°  <  <p  <  135°  just  before  all  of  the  pressures  increase  to  Cp  =  —  1  by  <p  =  135°  in 
the  steady  M=  0.4  case.  For  A cp  =  -28.2°,  Mach  number  is  approaching  the  minimum  value  as  a 
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increases.  Thus,  Cpxr  simultaneously  approaches  the  peak  values  for  the  cycle  (CPiCr  ~  -5.9).  A 
supercritical  condition  does  not  occur  for  this  phasing.  Greater  suction  is  developed  near  the 
leading  edge  until  cp  =  135°.  The  pressure  rise  which  indicates  stall  starts  with  the  upstream-most 
tap  and  proceeds  from  the  leading  edge  and  transpires  most  abruptly  for  the  A<p  =  -28.2°  case. 
The  phasing  alignment  in  the  second  case  (Atp  =  +46.4°)  is  such  that  Cpxr  =  -3.55  when  stall 
begins.  This  critical  pressure  is  approximately  equal  to  that  of  the  steady  M  =  0.4  case.  Recall 
that  the  Acp  =  +46.4°  and  steady  M  =  0.4  data  exhibited  moment  stall  at  the  exact  same  phase  in 
the  pitch  oscillation  cycle.  Three  upper  surface  taps  recorded  supercritical  pressures  in  the 
second  coupled-oscillation  case  (A (p  =  +46.4°)  and  the  ensuing  stall  was  relatively  gradual.  In  the 
final  coupled  oscillation  case,  for  A (p  =  +145.4°,  the  chordwise  extent  of  the  supercritical  region 
spans  four  taps  as  the  freestream  Mach  number  peaks  (Cpxr  =  -2.7)  and  the  airfoil  stalls. 
Comparing  the  rate  of  the  pressure  rise  between  all  four  cases  it  is  evident  that  as  the  extent  of 
the  supercritical  flow  on  the  airfoil  increases,  the  stall  process  becomes  increasingly  gradual  in 
nature. 
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A(|)  =  -28.2° - At|)  =  +46.4°  A(>  =  +145.4° - steady  M  =  0.4 

Figure  42:  Comparing  phase-averaged  upper  surface  (leading  edge)  Cp  with  CPlCr  for 
coupled-oscillation  and  steady  M  =  0.4  cases  and  k  =  0.0250. 

Only  two  cases  are  presented  in  Figure  43  (Cl  and  Cm)  and  Figure  44  (Cp  contours)  for  k  = 
0.0500.  The  phasing  of  the  coupled-oscillation  cases  is  a  close  match  to  the  target  with  a  leading 
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M  by  A<p  =  +13.3°.  This  case  is  compared  with  that  of  the  airfoil  executing  an  identical  pitch 
oscillation  profile  in  a  steady  M  =  0.4  airstream.  Apart  from  the  interval  for  90°  <  cp  <  180°,  the 
impact  of  time-varying  compressibility  on  Cl  and  Cm  hysteresis  is  marginal.  The  post-stall  and 
reattachment  process  during  the  downstroke  are  unaltered  even  for  this  highest  reduced 
frequency.  For  90°  <  (p  <  180°,  a  is  increasing  from  7.5°  to  22°.  The  lift  coefficient  and  lift  curve 
slope  are  higher  in  the  coupled-oscillation  case  throughout  this  interval.  The  difference  in  Cl 
between  the  two  cases  (A Cl)  increases  until  both  cases  reach  CiJnax  at  a  =  18.5°  with  a  maximum 
A  Cl  =  0.3.  The  slope  of  the  lift  stall  characteristic  is  approximately  equal  in  both  cases  indicating 
that  the  stall  type  is  unchanged.  Based  on  the  Cl  hysteresis  behavior,  the  manner  in  which  the 
dynamic  stall  process  is  modified  by  time-varying  compressibility  is  similar  for  all  three  reduced 
frequencies  for  which  the  phasing  is  similar  to  that  presented  in  Figure  37.  However,  for  k  = 
0.0125  and  0.0250,  coupled-oscillations  with  this  type  of  phasing  (A (p  ~  0°)  did  not  produce  a 
marked  differences  in  the  Cm  hysteresis.  The  peak  negative  Cm  increases  by  20%  for  the  k  = 
0.0500  cases  with  respect  to  the  steady  M=  0.4  result. 

Figure  44  juxtaposes  the  upper  surface  Cp  contours  for  the  two  k  =  0.0500  cases.  Suction 
development  preceding  the  LEV  fonnation  is  enhanced  for  the  coupled-oscillation  case.  The 
peak  is  Cp  =  -5.15  at  <p  =  139°  and  a  =  17°  for  the  coupled  oscillation  case,  while  the  peak  is  Cp 
=  -4  at  (p=  128°  and  a  =  15°  for  the  steady  M=  0.4  case.  As  with  cases  of  similar  phasing  for  k  = 
0.0125  and  0.0250,  the  magnitude  of  the  Cp  depression  created  by  the  LEV  signature  is 
amplified.  Observing  how  the  LEV  signature  is  modified  for  cases  of  similar  A<p,  it  seems  that 
the  LEV  amplification  relative  to  the  steady  M  =  0.4  cases  is  proportional  to  the  reduced 
frequency.  Along  with  the  exaggerated  low  pressure  signal,  the  chordwise  extent  of  the  LEV 
signature  relative  to  the  steady  M  =  0.4  data  indicates  a  stronger  LEV  as  it  convects  along  the 
upper  surface.  These  plots  do  not  necessarily  speak  to  the  rate  at  which  the  LEV  strength  decays, 
but  the  LEV  is  certainly  stronger  as  it  travels  downstream  past  the  quarter  chord,  thereby 
exacerbating  the  peak  negative  pitching  moment  as  shown  in  Figure  43. 
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Figure  43:  Cl  and  Cm  orbits  for  a  =  8.5°  ±  13°,  M  =  0.4  ±  0.08,  k  =  0.0500,  and 
representative  phasing  between  synchronous  oscillations. 
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Figure  44:  Cp  contours  for  a  =  8.5°  ±  13°,  M  =  0.4  ±  0.08,  k  =  0.0500,  and  representative 

phasing  between  synchronous  oscillations. 

The  results  presented  for  these  three  reduced  frequencies  suggests  that  compressibility  has  a 
first-order  effect  on  the  peak  Cl  and  the  onset  of  dynamic  stall,  while  the  acceleration  parameter 
plays  a  secondary  role  which  can  produce  lift  enhancement  especially  during  the  deceleration 
interval.  When  freestream  deceleration  coincides  with  the  upstroke  and  the  onset  of  dynamic 
stall,  it  precipitates  a  stronger  dynamic  stall  process  wherein  greater  suction  develops  on  the 
upper  surface,  Cimax  is  greater,  and  the  LEV  strength  is  amplified  along  with  the  pitching 
moment  spike.  The  hysteresis  effects  are  exaggerated  as  reduced  frequency  increases. 


J.  W.  Gregory,  W91 1NF1 1 10503  Final  Report,  p.  72 


4.  Conclusion 

An  advanced  dynamic  stall  wind  tunnel  experiment  was  conducted  in  which  the  free  stream 
velocity  was  made  to  oscillate  at  the  same  frequency  as  the  airfoil  pitch.  Pressure  contours,  along 
with  integrated  lift  and  moment  data,  were  presented  for  an  airfoil  subjected  to  various 
compressible  dynamic  stall  conditions.  The  dynamic  stall  problem  was  segmented  into  three 
parts  to  facilitate  understanding:  pitch  oscillations  in  a  fixed  compressible  freestream,  Mach 
oscillations  with  fixed  angle  of  attack,  and  coupled  pitch  and  Mach  oscillations. 

As  is  commonly  observed  with  steady  compressibility  effects,  the  max  Cl  in  dynamic  stall  is 
lowered  with  increased  freestream  M.  The  Cm  is  less  negative,  but  the  dimensional  moment 
magnitudes  are  obviously  still  much  higher,  and  fluctuations  are  very  significant.  One  critical 
aspect  of  compressibility  effects  is  a  change  from  leading-edge  dominated  dynamic  stall,  where 
the  inception  of  the  LEV  coincides  with  the  collapse  of  the  suction  peak,  to  trailing-edge 
dominated.  Suction  levels  near  the  leading  edge  are  sustained  as  the  freestream  Mach  number 
increases  whereas  the  suction  peak  completely  collapses  in  the  case  of  low  freestream  Mach 
number.  The  transfonnation  in  the  stall  process  from  leading-edge  to  trailing-edge  type  can 
inhibit  the  fonnation  and  strength  of  the  LEV  and  result  in  lower  peak  negative  moment 
magnitude. 

For  fixed  alpha,  oscillating  Mach  work,  the  steady-state  stall  angle  of  attack  for  the  mean 
Mach  number  of  the  oscillations  is  a  critical  parameter.  Both  lift  and  moment  tend  to  oscillate 
about  their  steady-state  flow  values  and  both  are  greater  than  steady-state  while  the  freestream  is 
decelerating  and  smaller  while  accelerating.  Reduced  frequency  exacerbates  this  hysteresis. 
Higher-frequency  periodic  fluctuations  appear  when  the  static  stall  angle  is  reached  (after  abrupt 
loss  of  leading-edge  suction  peak). 

The  phase  relationship  between  the  freestream  velocity  and  the  airfoil  angle  of  attack  was 
varied  and  found  to  have  significant  effects  on  the  aerodynamic  loads.  In  particular: 

1)  In-phase  oscillation  of  the  airfoil  pitch  and  the  free  stream  velocity  was  found  to  cause 
increased  lift  curve  slope  and  increased  stall  angle  compared  to  an  airfoil  with  oscillating  pitch  in 
steady  flow.  The  increased  stall  angle  is  associated  with  a  more  sudden  loss  of  lift. 

2)  Out-of-phase  oscillation  had  the  opposite  effect,  decreasing  both  the  lift  curve  slope  and 
the  stall  angle.  This  was  associated  with  a  more  gradual  stall  event  and  loss  of  lift. 
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3)  The  magnitude  of  peak  negative  moment  coefficient  was  found  to  be  stronger  for  in-phase 
motion  than  for  out-of-phase  motion,  as  the  leading  edge  vortex  convects  further  along  the  upper 
surface  of  the  model. 

4)  The  time  history  of  the  moment  coefficient  differs  mainly  quantitatively  between  phase 
angles,  but  retains  the  same  qualitative  shape;  this  contrasts  with  visible  qualitative  changes  to 
moment  orbits  at  different  static  Mach  numbers. 

Compressibility  has  a  first-order  effect  on  the  peak  Cl  and  the  onset  of  dynamic  stall,  while 
the  acceleration  parameter  plays  a  secondary  role  which  can  produce  lift  enhancement  especially 
during  the  deceleration  interval.  When  freestream  deceleration  coincides  with  the  upstroke  and 
the  onset  of  dynamic  stall,  it  precipitates  a  stronger  dynamic  stall  process  wherein  greater  suction 
develops  on  the  upper  surface,  CLmax  is  greater,  and  the  LEV  strength  is  amplified  along  with  the 
pitching  moment  spike.  The  hysteresis  effects  are  exaggerated  as  reduced  frequency  increases. 

Oscillation  of  freestream  velocity  is  a  fundamental  aspect  of  the  operating  environment  of  the 
lifting  surfaces  of  a  rotorcraft,  and  its  absence  in  wind  tunnel  studies  was  shown  to  result  in 
changes  in  the  measured  perfonnance  of  the  test  article.  Accurately  simulating  oscillations  of 
pitch  and  freestream  velocity  in  the  airfoil  environment  is  necessary  to  fully  reproduce  the 
aerodynamic  loads  experienced  by  a  rotor  in  forward  flight. 
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Appendix  1:  Airfoil  Coordinates 


Surface  Coordinates 

Surface  Derivatives 

Upper 

Lower 

Upper 

Lower 

x/c 

y/c 

y/c 

x/c 

dy/dx 

dy/dx 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00020 

0.00200 

-0.00145 

0.00010 

10.05190 

-7.30450 

0.00080 

0.00395 

-0.00287 

0.00050 

3.24740 

-2.36310 

0.00199 

0.00648 

-0.00457 

0.00140 

2.12040 

-1.42420 

0.00299 

0.00803 

-0.00545 

0.00249 

1.55160 

-0.87640 

0.00449 

0.00987 

-0.00645 

0.00374 

1.23030 

-0.66740 

0.00698 

0.01239 

-0.00770 

0.00573 

1.01220 

-0.50470 

0.00997 

0.01492 

-0.00888 

0.00848 

0.84550 

-0.39280 

0.01595 

0.01908 

-0.01070 

0.01296 

0.69460 

-0.30540 

0.02193 

0.02250 

-0.01218 

0.01894 

0.57240 

-0.24580 

0.02792 

0.02545 

-0.01345 

0.02493 

0.49240 

-0.21260 

0.03390 

0.02804 

-0.01459 

0.03091 

0.43370 

-0.19070 

0.03988 

0.03037 

-0.01563 

0.03689 

0.38940 

-0.17440 

0.04586 

0.03249 

-0.01659 

0.04287 

0.35520 

-0.16100 

0.05185 

0.03445 

-0.01749 

0.04885 

0.32680 

-0.14930 

0.05783 

0.03625 

-0.01831 

0.05484 

0.30090 

-0.13840 

0.06780 

0.03890 

-0.01957 

0.06281 

0.26620 

-0.12580 

0.07777 

0.04114 

-0.02069 

0.07278 

0.22470 

-0.11260 

0.08774 

0.04302 

-0.02171 

0.08275 

0.18790 

-0.10180 

0.09771 

0.04458 

-0.02264 

0.09272 

0.15720 

-0.09350 

0.11266 

0.04650 

-0.02391 

0.10519 

0.12850 

-0.08510 

0.12762 

0.04805 

-0.02506 

0.12014 

0.10360 

-0.07710 

0.14257 

0.04935 

-0.02612 

0.13510 

0.08630 

-0.07090 

0.15753 

0.05044 

-0.02710 

0.15005 

0.07340 

-0.06550 

0.17249 

0.05139 

-0.02801 

0.16501 

0.06290 

-0.06080 

0.18744 

0.05218 

-0.02885 

0.17996 

0.05350 

-0.05610 

0.20240 

0.05286 

-0.02963 

0.19492 

0.04510 

-0.05190 

0.21735 

0.05343 

-0.03034 

0.20987 

0.03800 

-0.04750 

0.23231 

0.05391 

-0.03099 

0.22483 

0.03230 

-0.04340 

0.24726 

0.05432 

-0.03158 

0.23978 

0.02750 

-0.03950 

0.27717 

0.05496 

-0.03259 

0.26222 

0.02130 

-0.03390 

0.30708 

0.05537 

-0.03340 

0.29213 

0.01370 

-0.02700 

0.33699 

0.05556 

-0.03401 

0.32204 

0.00650 

-0.02030 

0.37687 

0.05549 

-0.03451 

0.35693 

-0.00180 

-0.01250 

0.41675 

0.05504 

-0.03464 

0.39681 

-0.01140 

-0.00330 
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0.43669 

0.05466 

-0.03456 

0.42672 

-0.01890 

0.00400 

0.45664 

0.05418 

-0.03438 

0.44667 

-0.02410 

0.00920 

0.47658 

0.05360 

-0.03409 

0.46661 

-0.02940 

0.01450 

0.49652 

0.05290 

-0.03368 

0.48655 

-0.03490 

0.02020 

0.51646 

0.05209 

-0.03317 

0.50649 

-0.04040 

0.02600 

0.53640 

0.05118 

-0.03253 

0.52643 

-0.04600 

0.03170 

0.55634 

0.05015 

-0.03179 

0.54637 

-0.05150 

0.03720 

0.57628 

0.04901 

-0.03095 

0.56631 

-0.05720 

0.04220 

0.59622 

0.04776 

-0.03002 

0.58625 

-0.06290 

0.04670 

0.61616 

0.04638 

-0.02900 

0.60619 

-0.06890 

0.05090 

0.63610 

0.04488 

-0.02790 

0.62613 

-0.07550 

0.05510 

0.65604 

0.04322 

-0.02672 

0.64607 

-0.08300 

0.05940 

0.67598 

0.04139 

-0.02545 

0.66601 

-0.09170 

0.06380 

0.69592 

0.03937 

-0.02409 

0.68595 

-0.10150 

0.06820 

0.71586 

0.03714 

-0.02264 

0.70589 

-0.11170 

0.07250 

0.73580 

0.03472 

-0.02112 

0.72583 

-0.12140 

0.07630 

0.75574 

0.03214 

-0.01954 

0.74577 

-0.12940 

0.07930 

0.77568 

0.02945 

-0.01792 

0.76571 

-0.13510 

0.08140 

0.79562 

0.02668 

-0.01627 

0.78565 

-0.13860 

0.08260 

0.81556 

0.02387 

-0.01462 

0.80559 

-0.14090 

0.08300 

0.83550 

0.02101 

-0.01296 

0.82553 

-0.14340 

0.08320 

0.85544 

0.01809 

-0.01129 

0.84547 

-0.14660 

0.08370 

0.87538 

0.01509 

-0.00960 

0.86541 

-0.15020 

0.08480 

0.89532 

0.01205 

-0.00786 

0.88535 

-0.15250 

0.08700 

0.91526 

0.00905 

-0.00608 

0.90529 

-0.15070 

0.08940 

0.93520 

0.00623 

-0.00429 

0.92523 

-0.14130 

0.08990 

0.95514 

0.00385 

-0.00261 

0.94517 

-0.11930 

0.08420 

0.97508 

0.00229 

-0.00133 

0.96511 

-0.07830 

0.06440 

0.98506 

0.00199 

-0.00099 

0.98007 

-0.03020 

0.03340 

0.99503 

0.00214 

-0.00086 

0.99004 

0.01480 

0.01300 

1.00000 

0.00241 

-0.00080 

0.99751 

0.05480 

0.01210 
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Appendix  2:  Pressure  Tap  Locations 


Upper  Surface 

Tap  No. 

x/c 

y/c 

z/c 

1 

0.000282 

-0.00176 

0.172067 

2 

0.022804 

0.022955 

0.164312 

3 

0.039563 

0.03025 

0.158541 

4 

0.056403 

0.035838 

0.152743 

5 

0.073554 

0.040244 

0.146837 

6 

0.109371 

0.046116 

0.134504 

7 

0.129914 

0.048267 

0.127431 

8 

0.159442 

0.050573 

0.117264 

9 

0.180233 

0.051815 

0.110105 

10 

0.201131 

0.052807 

0.102909 

11 

0.222392 

0.053599 

0.095588 

12 

0.24413 

0.054243 

0.088103 

13 

0.266187 

0.054752 

0.080508 

14 

0.288719 

0.055143 

0.07275 

15 

0.311568 

0.05541 

0.064882 

16 

0.334733 

0.055562 

0.056906 

17 

0.368783 

0.055517 

0.045182 

18 

0.391475 

0.055387 

0.037368 

19 

0.414166 

0.055077 

0.029555 

20 

0.436855 

0.05466 

0.021743 

21 

0.459541 

0.054104 

0.013931 

22 

0.482224 

0.053409 

0.006121 

23 

0.510572 

0.052342 

-0.00364 

24 

0.539698 

0.051014 

-0.01367 

25 

0.571175 

0.049311 

-0.02451 

26 

0.604995 

0.047166 

-0.03615 

27 

0.64115 

0.044471 

-0.0486 

28 

0.679622 

0.041037 

-0.06185 

29 

0.72038 

0.036608 

-0.07588 

30 

0.764967 

0.030904 

-0.09124 

31 

0.815746 

0.023845 

-0.10872 
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Lower  Surface 

Tap  No. 

x/c 

y/c 

z/c 

1 

0.000282 

-0.00176 

0.172067 

32 

0.037661 

-0.01525 

0.159196 

33 

0.054657 

-0.01788 

0.153344 

34 

0.071558 

-0.02001 

0.147524 

35 

0.088494 

-0.02178 

0.141693 

36 

0.10608 

-0.02337 

0.135637 

37 

0.124466 

-0.02483 

0.129307 

38 

0.156017 

-0.02701 

0.118443 

39 

0.18906 

-0.02894 

0.107065 

40 

0.223693 

-0.03062 

0.09514 

41 

0.259913 

-0.03203 

0.082669 

42 

0.297718 

-0.03317 

0.069651 

43 

0.338683 

-0.03403 

0.055546 

44 

0.373208 

-0.03448 

0.043658 

45 

0.406301 

-0.03463 

0.032263 

46 

0.439394 

-0.03454 

0.020868 

47 

0.474061 

-0.03413 

0.008932 

48 

0.514371 

-0.03322 

-0.00495 

49 

0.561622 

-0.03158 

-0.02122 

50 

0.61878 

-0.02886 

-0.0409 

51 

0.686976 

-0.02471 

-0.06438 

52 

0.743046 

-0.02055 

-0.08369 

53 

0.797305 

-0.01613 

-0.10237 
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